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ABSTRACT 

Core-collapse supernovae are dramatic explosions marking the catastrophic 
end of massive stars. The only means to get direct information about the super- 
nova engine is from observations of neutrinos emitted by the forming neutron star, 
and through gravitational waves which are produced when the hydrodynamic flow 
or the neutrino flux is not perfectly spherically symmetric. The multidimension- 
ality of the supernova engine, which breaks the sphericity of the central core such 
as convection, rotation, magnetic fields, and hydrodynamic instabilities of the su- 
pernova shock, is attracting great attention as the most important ingredient to 
understand the long-veiled explosion mechanism. Based on our recent work, we 
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summarize properties of gravitational waves, neutrinos, and explosive nucleosyn- 
thesis obtained in a series of our multidimensional hydrodynamic simulations and 
discuss how the mystery of the central engines can be unraveled by deciphering 
these multimessengers produced under the thick veils of massive stars. 

Subject headings: supernovae: collapse — gravitational waves — neutrinos - 
hydrodynamics 



Introduction 



The majority of stars more massive than ~8 M & end their lives as core-collapse super- 
novae. They have long attracted the attention of astrophysicists because they have many 
facets playing important roles in astrophysics. They are the mother of neutron stars as well 
as black holes; they play an important role for acceleration of cosmic rays; they influence 
galactic dynamics triggering further star formation; they are gigantic emitters of neutrinos 
and gravitational waves. They are also a major site for nucleosynthesis, so, naturally, any 
attempt to address human origins may need to begin with an understanding of core-collapse 
supernovae (CCSNe). 

Current estim ates of CCSN rates in our Galaxy predict one event every ~ 40 ± 10 year 



( lAndo et al.l 120051 ) . When a CCSN event occurs in our galactic center, copious numbers 
of neutrinos are produced, some of which may be detected on the earth. Such "supernova 
neutrinos" will carry valuable information from deep inside the core. In fact, the detection 
of neutrinos from SN1987A (albeit in the Large Magellanic Cloud) opened up the Neutrino 



■ — 

Hirata et al. 


1987; 


Bionta et al. 


1987) 



events from SN1987A, these events have been studied extensively (yielding ~ 500 papers) 
and have all owed us to have a con fidenc e that our ba s ic pic ture of the supernova physics is 



correct (e.g., ISato fc Suzuki! (J1987J), see iRaffeltl ( 120121 . |2002| ) for a recent review). Recently 



significant pr ogress has been made in the l arge water Cherenkov detectors such as S uper- 
Kamiokande (ITotsukal Il992l ) and IceCube (jHurtqvist & IceCube collaboration! l201ll ). and 
also in the liquid scintillator detector as KamLAND (jSuzukil Il999l ) . If a supernova occurs in 
our Ga lactic center (~ 10 kpc), about 10 5 v e events are estimated to be detectable by IceCube 



e.g. : 



IceCube Collaboration et al.l ( 120111 ) and references therein) . Those successful neutrino 
detections are important not only to study the supernova physics but also to unveil the 
natur e of neutrinos itself such as the neutrino oscillation parameters and the mass hierarchy 



(.e.g. 



Raffeltl ( 120121 ) for a recent review). 
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CCSNe are now about to start even another astronomy, Gravitati onal- Wave Astronom y. 

Currently long-baseline laser interferometers such as LIGO (USA, e .g-.lAbbott et al.l (|2005l)). 

VIRG O (Italyfl GEO600 (Germanyfl and TAMA300 (Japan, e.g. JAndo fc the TAMA Collaboration 



(2005)) are currently operational and preparing for the first observation (see, e.g.. lHough et al 



(120051 ) for a recent review), by which the prediction by Einstein's theory of General Relativ- 
ity (GR) can be confirmed. These instruments are being updated to their Advanced status, 
and may start takin g data, possibly detecting GWs for the first time, as soon as 2015 (see 
van der Sluyd (120111 ) for a recent review). In fact, Advanced LIGO/VIRGO, which is an up- 
grade of the initial LIGO and VIGRO, are expected to be completed by 2015 and will increase 

the o bservable detection volume by a factor of ~ 1000 (IHarry fc the LIGO Scientific Collaboration! 

2010). The Large-scale Cryogenic Gravitational- wave Telescope (LCGT. lKuroda &: LCGT Collaboration 
(120 10l )) in Japan was funded in late 2010, which is being built under the Kamioka mine and 
is expected to take its first data in 2016. At such a high level of precision, those GW detec- 
tors are sensitive to many different sources, including chirp, ring-down, and merger phase s 
of black-hol e and neutron star binaries (e.g., ISathyaprakash fc Schutzl (I200S ); iFaberl ( 120091 ); 
Duez ( 2010 )). neutron star nor mal mode oscilla tions (e.g., Andersson et ali ( 2010)), rotating 



neutro n sta r moun t ains (e.g.. iHorowitzl (120111 )). and CCSN explosions (e.g., iKotake et al. 



(120061 ): lOttl (l2009bl ); iFryer fc New! (1201 ll ) for recent reviews), on the final of which we focus 
in this article. 



According to the Einstein's theory of GR (e.g., iShapiro &: Teukolskyl (119831 )). no GWs 
can be emitted if gravitational collapse of the supernova core proceeds perfectly spherically 
symmetric. To produce GWs, the gravitational collapse should proceed aspherically and 
dynamically. Observational evidence gathered over the last few decade s has pointed towards 
CCSNe indeed being generally aspherical (see I Wang fc Wheelerl (120081 ) for a recent review). 
The most unequivocal example is SN1987A. To explain the light-curve, a large amount of 
mixi ng of Ni outwar d to the H-He int e rface a nd of H i nward into the He-core was required 
fe.g.. lWooslevl Jl988h : lilinnikov et al.l fcoooh : lutrobir] tooi ). and lWooslev & Hegerl (hoOTf ) 
for a review). Such mixing processes coming from the Rayleigh- Taylor instability at the 
interface with differ ent composit i ons af t er shock pass a ge ha ve been examined exten sively 



so far by 2D fe .g.. lArnett et a 



3D simulations (jHammer et al 



Jl989h : iMiiller et al.l ( jl99lh : iKifonidis et all (12003^ and 
2010^ 7 The asymmetry of iron and ni ckel lines in SN1987A 



was proposed to be explained, if the ex plosion occurs in a jet-l ike ( jNagataki et al.l 119971 ; 
Nagatakill2000l ) or in a unipolar manner (jUtrobin fc Chugaill201ll ). More directly, the HST 
images of SN1987A are showing that the expanding envelope is elliptical with the long 



1 http: / / www.ego-gw.it / 
2 http://geo600. aei.mpg.de/ 
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axis aligned with the rotation axis inferred from the ring (IWang et al.l ( 120021 ). see how- 



ever 



Kjaer et al.l (120101 ) for a recent counter argument). The aspect ratio and position an- 
gle of the symmetry axis are consistent with those predicted earlier from the observations 
of speckle and linear polarization . Wh a t is more, the li near p olarization became greater 
as time passed (e.g., IWang et al.l (120011 ); iLeonard et al.l (l2006al Jbl)). a fact which has been 



used to argue that the central en g ine of the exp l osion i s responsible for the n on-sphericity 
Kifonidis et all tooA hood ): IScheck et all fcood ): burrows et all J2007cl laR By per- 



le-g. 

forming a series of time-de pendent, non-LTE ( local t hermodyna mic equilibrium), rad i ation 



transport simulations (e.g., bessart fc Hillierl fcoid bood . l2005l lR). bessart fc Hillierl (boilh 



recently pointed out that asymmetry in the eje cta can explain the increase in the continuum 



polarization observed at the nebular phases (ILeonard et al.l l2006d ; IChornock et al.ll2010l ). 



HST in the visual 


(IFesen et al. 


2006i 


2011) 


(Aschenbach et al. 


19951: 


Hughes et al. 


2000 



201lh and by ROSAT, C handra, and XMM-Newton 



Miceli et al.ll2008l ) in X-rays also provide evi- 



dence that small- and large-scale inhomogeneities (and perhaps even fragmentation) are a 
common feature in supernova explosions. 

Advancing ability of the HST has enabled the direct obs ervation of the progenitors of 
nearby CCSNe from pre-explosion images (see ISmarttl (120091 ) for a review). Although ob- 
servational me asurements of t he pr ogenitor masses are currently not many and still highly 
uncertain (see Murphy et al.l (1201 ll ) for collective references therein), evidence has accumu- 
lated that S N type H-pla teau (II-P) comes predominantly from stars in the range about of 
8-16 M Q ( ISmarttl 120091 ) . A gen eric explosion energy of the SN II-P in the mass range is 
roughly on the order of 10 51 erg ( lUtrobin fc Chugaill201ll ). however a large diversity of the 
explosion signatures (i.e. explosion ene rgy and the s ynthesized Ni mass) have been so far 
observed from quite similar progenitors (ISmarttl |2009| ). For example, the inferred 56 Ni mass 
and the kinetic energies differ by a factor of five between SNe 2005cs and 2003gd, both of 
them are among the golden events in which enough information was obtained to give an 
accurate estimate on a color or spectral type of the progenitor and the initial mass. More 
massive stars are expected to lose much of their mass and explode as hydrogen-stripped SNe 
(Ib/c a nd lib). Among them, th e type Ic-BL SNe, which are associated with long gamma-ray 
bursts ( IWoosley fc Blooml 120061 ). all show much broader lines than SNe Ic. Due to t he large 



kineti c energies of 2 — 5 x 10 52 erg, they have been referred to as "hypernovae" (e.g., iModjaz 
( 120111 ) for a recent review). These events are likely to come from interacting binaries in 
which the primary exploding star has a mass lower than what is usually associated w i th evo - 



lution to the massive WR phase (e.g.. iPodsiadlowski et al.l ( 120041 . |2010| ); iFryer et al.l ( 120071 ); 



SCO 



Kasen et al.l |2006h : lKasen fc Wooslevl (|2009h for an alternative approach for the light-curve modeling. 
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Smith et al.l ( 1201 ll ). see however I Yoon fc Langerl (I2005a[); Yoon et al.l (120101 ) ) . In addition to 
the two branches mentioned above, iNomoto et al.l (120061 ) predicted yet another branch, in 
which the SNe II-P with higher progenitor mass result in fainter explosions. By connecting 
to candidate SN explosion mechanisms or to progenitor structures, it is indeed best if one 
could obtain a unified picture to understand the mentioned wide diversity which is not only 
dependent on the progenitor masses but also on the evolution scenarios (a single vs. binary 
evolution). But, this is an area for future study firstly because it is too computationally ex- 
pensive to perform a long-term simulation that follows multidimensional (multi-D) dynamics 
consistently from the onset of core-collapse, through explosion, up to the nebular phased, 
secondly because we are still inaccessible to mu lti-D stellar evolution mod e ls, in which multi- 



D mo deling has been a major undertaking (see lArnett fc Meakinl (j2011al ); iMeakin Sz Arnett 
(120071 ) for recent developments). 



From a theoretical point of view, clarifying what makes the dynamics of the supernova 
engine deviate from spherical symmetry is essential in understanding the GW emission mech- 
anism. Here it is worth mentioning that GWs are primary observables, which imprint a live 
information of the central engine, because they carry the information directly to us with- 
out being affected in propagating from the stellar center to the earth. On the other hand, 
SN neutrinos are exposed to a number of (external) environmental effects, including self- 
interaction that induces co llective neutrino flavor oscillations predominantly in the vicinity 



of the neutrino sphere (see lRaffelt fc Smirnovl 120071 ; iDuan et al.ll2008l . I201CH : iDasguptal 12010 



for reviews of the rapidly growing research fi eld and collective reference s therein) and the 
Mikheyev-Smirnov-Wolfenstein (MSW) effect (IMikheev fc Smirnovlll986l ) both in the stellar 
envelope and in the earth (see iRaffeltl (120121 ) for a review). Although the time profiles of 
neutr ino signals can b e pote ntially used like a tomography to monitor the envelope profile 
(e.g., IDuan fc Knellerl ( 120091 ) for a recent review), SN neutrinos generally could provide a 
rather indirect information about the central engine compared to GWs. 

The breaking of the sphericity in the supernova engine has been considered as the most 
important ingredient to understand the explosion mechanism, for which supernova theorists 
have been continuously keeping their efforts for the past ~ 40 years. Currently multi-D sim- 
ulations based on refined numerical models have shown several promising scenarios. Among 
the candidates are the neutrino heating me chanism aided by c onvection and hydrodynamic 
instabilities of the supernova sho ck (e.g., IJanka et al.l ( 120071 ) for a review), the acoustic 
mechanism (iBurrows et al.ll2007bl ). or the magnetohydrodynamic (MHD) mechanism (e.g., 
Kotake et alT ( 2006 ) see references therein). To pin down the true answer among the can- 



for example, see 
for binary models. 



Dessart et al. 



([201 ll ) for the most up-to-date ID modeling of spectra and light curves 
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didate mechanisms, GW signatures albeit being a primary observable, will not be solely 
enough and a careful analysis including neutrinos and photons should be indispensable. The 
current neutrino detectors is ready to broadcast the alert to astronomers t o let them know 



the a rrival of neutrinos (e.g., Supernova Early Warning Systems (SNEWS) (lAntonioli et al. 



2004 )). In addition to optical observation s using largest 8-10 m telescopes such as VLT and 
Subaru telescope (e.g., iPatat et al.l (120091 ); iTanaka et al.l ( 120091 )). the planned thirty- meter- 
telescope (TMT)@ with a refined spectropolarimetric technique is expected to detect more 
than 20 events of the SN polarization per yeai(§. This should provide valuable information 
to understand the SN asymmetry with increasing statistics. Not only for understanding the 
origin of non-spherical ejecta morphologjfj] but more importantly for understanding the ori- 
gin of heavy elements, it is of crucial importance to accurately determine nucleosynthesis in 
the SN ejecta, which naturally re quires a multi-D numerical modeling (see, IWoosley fc Heger 
( 2007 ); Thielemann et al. ( 2011 ) for recent reviews)0. 



Putting things together, the multidimensionality determines the explosion mechanism, 
in turn we may extract the information that traces the multidimensionality by the SN mul- 
timessengers, which would be only possible by a careful analysis on GWs, neutrinos, and 
photons. In this review article, we hope to bring together various of our published and unpub- 
lished findings from our recent multi-D supernova simulations and the obtained predictions of 
the SN multimessengers so far (for other high-en ergy astrophysical s ources such as magnetars, 



gamm a -ray bursts, and coa l escing binaries. seelAndo et al.l (120121); IChassande-Mottin et al. 



(2011 


); 


Marka et al. 


(2011 


); 



Aso et al 



( 120081 ) for 



recent reviews). Before we go into details from the next sections, we first have to draw a 
caution that the current generation of simulation results that we report in this article should 
depend on the next generation calculations by which more sophistication can be made not 
only in determining the efficiency of neutrino-matter coupling (the so-called neutrino trans- 
port calculation), but also in the treatment of general relativity. Therefore we provide here 
only a snapshot of the moving (long-run) documentary film whose headline we (boldly) chose 
to entitle as "multimessengers from CCSNe to bridge theory and observation" . 



Among the mentioned candidate mechanisms, we focus on the neutrino-heating mech- 



5 http://www. tmt.org/ 

6 M.Tanaka in private communication. 



Especially for nearby event 

8 It is worth mentio ning that radioactive de cay can affect the shape of the light curve for some peculiar 
SN 1987A-like events (jutrobin fc Chugai 2011 ). in which explosive nucleosynthesis plays an important role 
for the light-curve modeling. 



-7- 



anism in section [2] and the MHD mechanism in section [3X respectively. In each section, 
we first briefly summarize the properties of the explosion dynamics and then move on to 
discuss possible properties of the SN multimessengers paying particular attention to their 
detectability. It may be best if we can cover these SN messengers once for all in this arti- 
cle, but unfortunately not. What we have studied so far is limited to GWs and explosive 
nucleosynthesis in the neutrino-heating mechanism, and to GWs and neutrino signals in the 
MHD mechanism. To compensate the uncovered fields, the related references will be given. 
Although a number of excellent reviews already exist on various topics in this article, this 
one might go beyond such reviews by its new perspectives on the multimessenger astronomy. 



2. Neutrino-heating mechanism 



CCSN simulations have been counted as one of the most challenging subjects in compu- 
tational astrophysics. The four fundamental forces of nature are all at play; the collapsing 
iron core bounces due to strong interactions; weak interactions determine the energy and 
lepton number loss in the core via the transport of neutrinos; electromagnetic interactions de- 
termine the properties of the stellar gas; GR plays an important role due to the compactness 
of the proto-neutron star and also due to high velocities of the collapsing material outside. 
Naturally, such physical richness ranging from a microphysical scale (i.e. femto-meter scale) 
of strong/weak interactions to a macrophysical scale of stellar explosions has long attracted 
the interest of researchers, necessitating a world-wide, multi-disciplinary collaboration to 
clarify the theory of massive stellar core-collapse and the formation mechanisms of compact 
objects. 



Ever since the first num erical simulat i on of such events (IColgate fc White! Il966l ) , the 
neutrino-heating mechanism (Wilson 1985 ; Bethe fc Wilson 19851 ; Bethel 199oh . in which a 
stalled bounce shock could be revived via neutrino absorption on a timescale of several hun- 
dred milliseconds after bounce, has been the working hypothesis of supernova theorists for 
these ~ 45 years. However, the simplest, spherically-symmetric (I D) form of this mechanism 
fails to blow up canoni c al massive stars Rampp fc Janka ( 200ol ); Liebendorfer et al. ( 2001 ); 



Thompson et al.l (120031 ) ; ISumivoshi et al 



blast morphology (e.g., IWang fc Wheeler 



2005a). Pushed by mounting observations of the 
20081 ) mentioned above, it is now almost certain 



that the breaking of the spherical symmetry is the key to solve the supernova problem. So 
far a number of multi-D hydrodynamic simulations have been reported , which demo n strate d 
that hydrodyn a mic m o tions associated with c o nvective over t urn ( e. g.. iHerant et al.l ( 11994] ); 



Burrows et all (I1995f ): Ijanka fc Mullerl ( Il996h : iFrver et al.l (120021): iFrveri (l2004af) N ) as well 



as the Standing-Accretion-Shock-Instability (SASI, e.g.. iBlondin et al.l (120031 ); 



Scheck et al. 



(hooi l2006h: lOhnishi et al.l (bood . [2007h: iForiizzo et all ( ho06h : 



Fernandez fc Thompsonl (l2009bl Jal) ; llwakami et al.l (12008 



20091 ): 



Murphv fc Burrows! (120081 ) : 
Fernandez (j2010 ) see refer- 



ences therein) can help the onset of the neutrino-driven explosion. 

To test the neutrino heating mechanism in the multi-D context, it is of crucial impor- 
tance to solve accurately the neutrino-matter coupling in spatially non-uniform hydrody- 
namic environments. For the purpose, one ultimately needs to solve the six-dimensional 
(6D) neutrino radiation transport problem (three in space, three in the momentum space 
of neutrinos), which is a main reason why the supernova simulations stand out from other 
astrophysical simulations due to their complexity. In the final sentence of the last para- 
graph, we wrote "demonstrated" because the neutrino heating was gi ven by han d as a n 
input parame ter in most of the simulations cited above (see, however iFryer et al.l ( 120021 ); 
Fryerl (j2004al )). The neutrino heating proceeds dominantly via the charged current interac- 
tions (u e + n <=^ e~ + p, z/ e + p e + + n) in the gain region. The neutrino heating rate 
in the gain region can be roughly expressed as Q + oc L u (fi u )~ 1 /r 2 where L v is the neutrino 
luminosity emitted from the surface of neutrino sphere and it determines the amplitude of 
the neutrino heating as well as cooling , and r and ( fi„) is the distance from the stellar center 
and the flux factocl, respectively (e.g.. Ijankal ( 1200 ll )). For example, L v i s treat ed as an input 



parameter in the so-called 'light-bulb" approach (e.g., Janka fc Muller ( 1996h). This i 



of the most prevailing approximations in recent 3D simulations ( llwakami et al.ll2008l . 12009 



is one 



Wongwathanarat et al.l |2010| ; iNordhaus et al.l 120101 ) because it is handy to study multi-D 
effects on the neutrino heating mechanism (albeit on the qualitative grounds). To go beyond 
the light-bulb scheme, L v should be determined in a self-consistent manner. For the purpose, 
one needs to tackle with neutrino transport problem, only by which energy as well as angle de- 
pendence of the neutrino distribution function can be determined without any assumptions. 
Since the focus of this review is on the SN multi-messengers, a detailed discussion of various 
approximations and numerical techniques taken in the recent radiation-hydrodynamic SN 
simulations cannot be provided. Table 1 is not intended as a comprehensive compilation, 
but we just want to summarize milestones that have recently reported the neutrino-driven 
exploding models so far. 

In Table 1, the first column ("Progenitor") shows the progenitor model employed in each 
simulation. The abbreviation of "N H", "WHW", and "WW " means iNomoto h Mashimoto 



(ll988l ). IWooslev et al.l (120021 ) . and lWooslev & Weaver! (Il995l ). The second column shows SN 
groups with the published or submitted year of the corresponding work. "MPA" stands for 
the CCSN group in the Max Planck Institute for Astrophysics led by H.T. Janka and E. 



9 This quantity represents the degree of anisotropy in neutrino emission; (/i„) ~ 0.25 near at the neutrino 
sphere, = 1 in the free-streaming limit (r — ¥ oo). 
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Miiller. "Princeton+" stands for the group chiefly consisting of the staffs in the Prince- 
ton University (A. Burrows, J. Murphy), Caltech (C.D.Ott), Hebrew University (E. Livne), 
Universite de Provence (L.Dessart), and their collaborators. The SN group in the Basel 
university is led by M. Liebendorfer and F.K. Thielemann. " OakRidge+" stands for the SN 
group mainly consisting of the Florida Atlantic University (S. Bruenn) and the Oak Ridge 
National Laboratory (A. Mezzacappa, O.E.B. Messer) and their collaborators. Tokyo+ is 
the SN group chiefly consisting of the staffs in the National Astronomical Observatory of 
Japan (myself, T. Takiwaki), Kyoto university (Y. Suwa), Numazu College (K. Sumiyoshi), 
Waseda University (S. Yamada), and their collaborators. The third column represents the 
mechanism of explosions which are basically categorized into two (to date), namely by 
the neutrino-heating mechanism (indicated by " //-driven") or by the acoustic mechanism 
("Acoustic"). "Dim." in the fourth column is the fluid space dimensions which is one-, two-, 
or three-dimension (1,2, 3D). The abbreviation "N" stands for 'Newtonian,' while "PN"— 
for 'Post-Newtonian' — stands for some attempt at inclusion of general relativistic effects, 
and "GR" denotes full relativity. t exp in the fifth column indicates an approximate typical 
timescale when the explosion initiates and E exp represents the explosion energy normalized 
by Bethe(=10 51 erg) given at the postbounce time of t ph , both of which are attempted to be 
sought in literatureg^j. In the final column of " v transport" , " Dim" represents v momentum 
dimensions and the treatment of the velocity dependent term in the transport equations is 

symbolized by 0(v/c). The definition of the "RBR", "IDSA", and "MGFLD" will be given 
o n D I iie.ito ^hg^^ga limit of this article, we have to start th e story only after 2006 (see, e.g. 



Janka et al.l ( 120071 ) for a complete review, and also ICardalll ( 120051 ) for a similar table before 



2006). The first news of the exploding model was reported by the MPA group. By perform- 
ing radiation-neutrino-hydrodynamic simulations which includes one of the best available 
neutrino transfer approximations, they reported ID and 2D explosions for the 8.8 M Q star 
( jKitaura et al.ll2006l ; I Janka et al.l 120081 ) whose progenitor has a very tenuous outer envelope 
with steep density gradient (a characteristic property of AGB stars). Also in ID, the Basel-I- 
group reported explosions for 10 and 15 M & progenitors of WHW02 triggere d by the hy- 



pothe sised first-order QCD phase transition in the protoneutron star (PNS) ( ISagert et al. 



20091 ). To date, these two are the only modern numerical results where the neutrino-driven 



mechanism succeeded in ID. In the 2D MP A simulations, they obtained explosions for a 
non-rot ating 11.2 progenit or of WHW02 ( iBuras et al.ll2006ai ). and then for a 15M pro- 
genitor ( Marek fc Janka 2009 ) of WW95 with a relatively rapid rotation imposed"! They 
newly brought in the so-called "ray- by-ray" approach (indicated by "RBR" in the table), in 



10 



but if wc cannot find them, we remain them as blank 



1 by comparing the pre c ollapse angular velo city to the one predicted in a recent stellar evolution calculation 



Maeder fc Mevnetl (|2000l l: iHeger et al.l (|2005() . 
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Progenitor 


Group 


Mechanism 


Dim. 


^exp 


-^exp(B) 


v transport 




(Year) 




(Hydro) 


(ms) 


@t ph (ms) 


(Dim, 0(v/c)) 




MPA 


//-driven 


1D(2D) 


-200 


0.1 


Boltzmann 


8.8 M 


(2006,2011) 




(PN) 




(~800) 


2, 0(v/c) 


(NH88) 


Princeton+ 


//-driven 


2D 


<125 


0.1 


MGFLD 




(2006) 




(N) 




- 


1, (N) 


10 M 


Basel 


z/+(QCD 


ID 


255 


0.44 


Boltzmann 


(WHW02) 


(2009) 


transition) 


(GR) 




(350) 


2, (GR) 


11 Mq 


Princeton+ 


Acoustic 


2D 


>550 


-0.1* 


MGFLD 


(WW95) 


(2006) 




(N) 




(1000) 


1, (N) 




MPA 


//-driven 


2D 


-100 


~ 0.005,0.025 


"RBR" Boltz- 


11.2 Mq 


(2006,2012) 




(PN,C-GR) 


-£00 


-200,000 


mann, 2, 0(v/c) 


(WHW02) 


Princeton+ 


Acoustic 


2D 


>1100 


-0.1* 


MGFLD 




(2007) 




(N) 




(1000) 


1, (N) 




Tokyo+ 


i/-driven 


3D 


-100 


0.01 


IDSA 




(2011) 




(N) 




(300) 


1, (N) 


12 Mq 


Oak Ridge + 


i/-driven 


2D 


-300 


0.3 


"RBR" MGFLD 


(WHW02) 


(2009) 




(PN) 




(1000) 


1, 0(v/c) 


13 Mq 


Princeton+ 


Acoustic 


2D 




-0.3* 


MGFLD 


(WHW02) 


(2007) 




(N) 




(1400) 


1, (N) 


(NH88) 


Tokyo+ 


//-driven 


2D 


-200 


0.1 


IDSA 




(2010) 




(N) 




(500) 


1, (N) 


15 Mq 


MPA 


//-driven 


2D 


-600 


0.025,0.125 


Boltzmann 


(WW95) 


(2009,2012) 




(PN,C-GR) 


~400 


(-700,500) 


2,0(v/c) 


(WHW02) 


Princeton+ 


Acoustic 


2D 


- 


- 


MGFLD 




(2007) 




(N) 




(-) 


1, (N) 




OakRidge+ 


i/-driven 


2D 


-300 


- 0.3 


"RBR" MGFLD 




(2009) 




(PN) 




(600) 


l,0(v/c) 


20 Mq 


Princeton+ 


Acoustic 


2D 


>1200 


-0.7* 


MGFLD 


(WHW02) 


(2007) 




(N) 




(1400) 


1, (N) 


25 Mq 


Princeton+ 


Acoustic 


2D 


>1200 




MGFLD 


(WHW02) 


(2007) 




(N) 




(-) 


1, (N) 




Oak Ridge+ 


i/-driven 


2D 


-300 


- 0.7 


"RBR" MGFLD 




(2009) 




(PN) 




(1200) 


1, 0(v/c) 



Table 1: Selected lists of recent milestones reported by SN groups around the world 
("Group"), which obtained explosions by the neutrino- heating mechanism (indicated by 
"i/-driven") or the acoustic mechanism ("Acoustic") (See text for more details). 
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which the neutrino transport is solved along a given radial direction assuming that the hydro- 
dynamic medium for the direction is spherically symmetric. This method, which reduces the 
2D p roblem partly to ID, fits well with their original ID Boltzmann solver (IRampp &: Janka 



2000 



12 



For 2D hydrodynamic simulations with the ray- by-ray transport, one needs to solve 
the 4D radiation transport problem (two in space and two in the neutrino momentum space). 
Regarding the explosion energies obtained in the MPA simulations, their values at their final 
simulation time are typically underpowered by one or two orders of magnitudes to explain 
the canonical supernova kinetic energy (~ 10 51 erg). But the explosion energies presented 
in their figures are still growing with time, and they c ould be as h i gh as 1 B if they were 
able to follow a muc h longer evolution as discussed in iBuras et al.l (j2006al ). Very recently, 
Mueller et al.l ( 120121 ) reported that more energetic explosions are obtained for the 11.2 and 
15M stars i n their GR 2D simula tions compared to the corresponding post-Newtonian 



models (e.g., Marek Sz Janka ( 20091 )). In the table, the data in italic character represent 



their most up-to-dat e results based on the 2D G R simulations using conforma lly-flatness 
approximation (e.g., iDimmelmeier et al.l (J2002bJ); ICordero- Carrion et al.l (120091 ). indicated 
by " C-GR n in the table). 

It is rather only recently that fully 2D multi-angle Bolt z mann transport sim ulations 
become practicable by the Princeton+ group (lOtt et al.ll2008l ; iBrandt et al.l 1201 ll ). In this 
case, one needs to handle the 5D problem for 2D simulations (two in space, and three in 
the neutrino momentum space). However this scheme is very computationally expensive 
currently to perform long-term supernova simulations. In fact, the most recent 2D work by 
Brandt et al.l (120111 ) succeeded in following the dynamics until ~ 400 ms after bounce for a 
non-rotating and a rapidly rotating 20 M model of WHW02, but explosions seemingly have 
not been obtained in such an earlier phase either by the neutrino-heating or the acoustic 
mechanism. 

In the table, "MGFLD" stands for the Multi-Group Flux-Limited Diffusion scheme 
which e li mina tes the angular dependence of the neutrino distribution function (see, e.g., 
Bruennl (119851 ) for more details). For 2D simulations, one needs to solve the 3D prob- 
lem, namely two in space, and one in the neutrino momentum space. By implementing 



the M GFLD algorithm to the CHIMERA code in a ray-by-ray fashion (e.g., iBruenn et al. 
(120 10l )). Bruenn et al. (2009) obtained neutrino-driven explosions for non- rotating progen- 
itors in a relatively wide range in 12, 15, 20, 25 M of WHW02 (see table). These models 
tend to start exploding at around 300 m after bounce, and the explosion energy for the 



12 Note that the ray-by-ray approach has an advantage compared to other approx i mation schemes, such 
that it can fully take into account the available neutrino reactions (e.g., Buras et al. ( 2006bl) for references 
therein) and also give us the most accurate solution for a given angular direction. 
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l ongest running mo del of the 25 M progenitor is reaching to 1 B at 1.2 s after bounce 
feruenn et aDhoiof ). 



On the other ha nd, the 2D MGFLD simulations implemented in the VULCAN code 
( jBurrows et al.l 120061 ) obtained explosions for a variety of progenitors of 11, 11.2, 13,15, 20, 
and 25 M & not by the neutrino-heating mechanism but by the acoustic mechanism! 13 !. The 
acoustic mechanism relies on the revival of the stalled bounce shock by the energy depo- 
sition via the acoustic waves that the oscillating protoneutron stars (PNSs) would emit in 
a much delayed phase (~ 1 second) compared to the conventional neutrino-heating mech- 
anism (~ 300 — 600 millisec onds). If the core pulsation energy given in Burrows et al. 
(2007) (jBurrows et al.ll2007d ) could be used to measure the explosion energy in the acoustic 
mechanism, they reach to 1 B after 1000 ms after bounce. 

By performing 2D simulations in which the spectral neu trino transport was solve d by 
the isotropic diffusion source approximation (IDS A) scheme (ILiebendorfer et al.ll2009l ). the 
Tokyo+ group reported explosions for a non-rotating and rapidly rotating 13 M progenitor 
of NH88. They pointed out that a stronger explosion is obtained for the rotating model 
comparing to the corresponding non-rotating model. The IDSA scheme splits the neutrino 
distribution into two components (namely the streamin g and trapped neu t rinos) , both of 
which are solved using separate numerical techniques (see ILiebendorfer et al.l (120091 ) for more 
details). The approximation level of the IDSA scheme is basically the same as the one of the 
MGFLD. The main advantage of the IDSA scheme is that the fluxes in the transparent region 
can be determined by the non-local distribution of sources rather than the gradient of the 
local intensity like in MGFLD. A drawback in the current version of the IDSA scheme is that 

v T and their anti-particles) as well as the ener gy-coupling 



v v , i.e., z/, 



heavy lepton neutrinos 

weak interactions have yet to be implemented. Extending the 2D modules in (jSuwa et al. 



2010) to 3D, they recently report ed explosions in the 3D models for an 11.2 M progenitor 
of WHW02 (ITakiwaki et al.ll2012l ). By comparing the convective motions as well as neutrino 
luminosities and energies between their 2D and 3D models, they pointed out whether 3D 
effects would help explosions or not is sensitive to the employed numerical resolutions. They 
argued that next-generation supercomputers are at least needed to draw a robust conclusion 
of the 3D effects. 

Having summarized a status of the current supernova simulations, one might easily 
see a number of issues that remain to be clarified. First of all, the employed progenitors 
usually rather scatter (e.g., Table 1). Different SN groups seem to have a tendency to employ 
different progenitors, providing different results. By climbing over a wall which may have 



13 For the 8.8 Mq progenitor, they obtained neutrino-driven explosions (see table 1). 



13 



rather separated exchanges among the groups, a detailed comparison for a given progenitor 
needs to be done seriously in the multi-D results (as have been conducte d in the Boltzmann 
ID s imulations between the MPA, Basel+, and Oak Ridge+ groups (ILiebendorfer et al. 
2005h ). 



In the nex t section, we briefly s ummarize the findings obtained in our 2D (ISuwa et al. 



20101 ) and 3D (ITakiwaki et al.l 120121 ) simulations, paying particular attention to how mul- 
tidimensionality such as SASI, convection, and rotation could affect the neutrino-driven 
explosions. 




100 150 200 250 300 350 400 450 
Time after bounce [ms] 



Fig. 1. — Time evolution of ID (thin gray lines) or 2 D (thin orange lines) hydro dynamic 
simulation ( ISuwa et al.l 120101 ) of a 13 M Q progenitor (INomoto fc Mashimotolll988l ). Thick 
lines in red (for 2D) and black (for ID) show the position of shock waves, noting for 2D that 
the maximum (top) and average (bottom) shock position are shown. The red dashed line 
represents the position of the gain radius, which is similar to the ID case (not shown). 



2.1. Multidimensionality in multi-D radiation hydrodynamic simulations 

Figured] depicts the difference between the time evolutions of ID (thin gray lines) or 2D 
(thin orange lines) simulation of the 13 M progenitor model. Until ~ 100 ms after bounce, 
the shock position of the 2D model (thick red line) is similar to the ID model (thick black 
line). Later on, however, the shock for 2D does not recede as for ID, but gradually expands 
and reaches 1000 km at about 470 ms after bounce. Comparing the position of the gain radius 
(red dashed line) to the shock position for ID (thick black line) and 2D (thick red line), one 
can see that the advection time of the accreting material in the gain region can be longer in 
2D than in ID. This longer exposure of cool matter in the heating region to the irradiation 
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Fig. 2. — Snapshot of the distribution of entropy (left half) and the ratio of the advection to 
the heating tirnescale (right half) for models o f the ID (le f t) and 2D (right) models at 200 
ms after bounce. These figures are taken from ISuwa et al.l ( 2010 ). 



< 10 



g 10 




200 300 
Time after Bounce [ms] 



Fig. 3. — SASI activity for the rotating model versus postbounce time. Shown are the 
coefficients of the dipole (£ = 1) and quadrupole (£ = 2) modes of the spherical harmonics 
of the aspherical s hock position, normalized to the amplitude of the i = mode (taken from 
Suwa et all (buluh ). 
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Fig. 4. — Time evolution of the explosion energy versus postbounce time for 2D models 
with and without rotation (left panel). The explosion energy is defined as the total energy 
(internal plus kinetic plus gravitational), integrated over all matter where the sum of the 
corresponding specific energies is positive. Models with "-hr" indicates the ones with higher 
numerical resolution, in which the mesh numbers in the lateral direction are doubled. Middle 
and right panels are snapshots of the density (left half) and the entropy (right half) for ID 
(left) and 2D rotating (right) model at the epoch when the shock reaches 1000 km. These 
figures are taken from lSuwa et al.l (120101 ). 
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of hot outstreaming neutrinos from the PNS is essential for the increased efficiency of the 
neutrino heating in multi-D models. This can be also depicted in Figure [2j The right-half 
shows T adv /r heat , which is the ratio of the advection to the neutrino heating timescaldij. For 
the 2D model (right panel), it can be shown that the condition of T adv /r hcat > 1 is satisfied 
behind the aspherical shock, which is deformed predominantly by the SASI (£ = 2 mode at 
the snapshot), while the ratio is shown to be smaller than unity in the whole region behind 
the spherical standing accretion shock (left pa neblD). Here SASI that develops due t o the 



advection-acoustic cycle in the supernova core (iFoglizzo fc Taggerll2000l ; lFoglizzoll200ll ). is a 
uni- and bipolar sloshing of the stalled bounce shock with pulsational strong expansion and 
contraction (seen as oscillations in the red curves in Figure 1). Comparing the left-half of 
each panel, the entropy for the 2D model is shown to be larger than for the ID model. This 
is also the evidence that the neutrino heating works more efficiently in 2D. 

To see clearly the effect s of stellar rotati on on the postbounce evolution, a rapidly 



rotating model was taken in ISuwa et al.l (120101 ). in which a constant angular frequency of 



Qq = 2 rad/s is imposed inside the iron core with a dipolar cut off (oc r~ 2 ) outside. This 
corresponds to (3 ~ 0.18% with (3 being the ratio of the rotational to the gravitational 
energy when the simulation starts. For the rotating model, the dominant mode of the shock 
deformation after bounce is almost always the i = 2 mode for the rotating model (e.g., 
Figure [31 although the t = 1 mode can be as large as the 1 = 2 mode when the SASI enters 
the non-linear regime (> 200 ms after bounce). In contrast to this rotation-induced 1 = 2 
deformation, the £ = 1 mode tends to be larger than the t = 2 mode for the 2D models 
without rotation in the saturation phase. 

The left panel of Figure H] shows the comparison of the explosion energies for the 2D 
models with and without rotation (models with "-rot" or "-2D", respectively). Although the 
explosion energies depend on the numerical resolutions quantitatively, they show a continuous 
increase for the rotating models. The explosion energies for the models without rotation, 
on the other hand, peak around 180 ms when the neutrino-driven explosion sets in (see also 
Figure 1), and show a decrease later on. The reason for the greater explosion energy for 
models with rotation is due to the bigger mass of the exploding material. This is because 
the north-south symmetric (£ = 2) explosion can expel more material than for the unipolar 
explosion (compare the middle and right panels in Figure H]). The explosion energies when 
we terminated the simulation are less than < 10 50 erg for all the computed models. For the 
rotating models, we are tempted to speculate that the explosion energies could increase later 



14 This quantity is known as a useful quantity to diagnose the success (T a d v /Theat ^ b i- e -, the neutrino- 
heating timescale is shorter than the adve ction timescale of materia l in the gain reg i on) or failure (r a .H y / Th Ra ,t ^ 
1) of the neutrino-driven explosion (e.g.. iBurrows fc Goshyl (|1993l ): IJankal (|200ll ): iThompson et al.1 (J2005|)). 
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on by a linear extrapolation. However, in order to identify the robust feature of an explosion, 
a longer-term simulation with improved input physics is needed. In combin ation with the 
assume d rapid rotation, magnetic fields should be also taken into account as in lBurrows et al. 
( )2007at ). which we are going to study as a follow-up of our rotating 2D models (Suwa et al. 
in preparation). 

Extending our 2D modules ment i oned above, we are currently running 3D simulations 
for an 11.2 M Q star of IWoosley et all (120021 ) with spectral neutrino transport that is solved 
by the IDSA scheme in a ray- by-ray manner (ITakiwaki et all 120121 ) . We briefly summarize 
the results in the following. 



10000 F 



1000 - 
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Time after bounce [ms] 
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Time after bounce [ms] 
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Fig. 5. — Time evolution of our 3D model (ITakiwaki et al.ll2012l ). visualized by mass shell 
trajectories in thin gray lines (left panel). Thick red lines show the position of shock waves, 
noting that the maximum (top), average (middle), and the minimum (bottom) shock position 
are shown, respectively. The green line represents the shock position of the ID model. "1.30" 
and "1.40" indicates the mass in unit of M Q enclosed inside the mass-shell. Right panel 
shows the evolution of average shock radii for the 2D (green line), and 3D (red line) models. 
The "3D low" (pink line) corresponds to the low resolution 3D model, in which the mesh 
numbers are ta k en to be half of the standard model (see text). These figures are taken from 



Takiwaki et all (I2012h . 



The left panel of Figure [5] shows mass-shell trajectories for the 3D (red lines) and ID 
model (green line), respectively. At around 300 ms after bounce, the average shock radius 
for the 3D model exceeds 1000 km in radius. On the other hand, a n explosion is not obtained 
for the ID model, which is in agreement with lBuras et all ( l2006af ). The right panel of Figure 
|5] shows a comparison of the average shock radius vs. postbounce time. In the 2D model, 
the shock expands rather continuously after bounce. This trend is qualitatively consistent 
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with the 2D result by iBuras et all (j2006al ) (see their Figure 15 for model sll2_128_f]0 



Comparing the shock evolution between our 2D (green line in the right panel of Figure 
[5]) and 3D model (red line), the shock is shown to expand much faster for 2D. The pink 
line labeled by "3D low" is for the low resolution 3D model, in which the mesh numbers are 
taken to be half of the standard model. Note that the 3D computational grid consists of 300 
logarithmically spaced, radial zones to cover from the center up to 5000 km and 64 polar (6) 
and 32 azimuthal ((f)) uniform mesh points, which are used to cover the whole solid angle. 
The low resolution 3D model has one-half of the mesh numbers in the direction (n^=16), 
while fixing the mesh numbers in other directions. Comparing with our standard 3D model 
(red line), the shock expansion becomes less energetic for the low resolution model (later 
than ~ 150 ms). Above results indicate that explosions are easiest to obtain in 2D, followed 
in order by 3D, and 3D (low). At first sight, this may look contradicted with the finding by 



Nordhaus et al.l (120101 ) who pointed out that explosions could be more easily obtained in 3D 



than in 2D. The reason of the discrepancy is summarized shortly as it follows. 

Figure [6] compares the blast morphology for our 3D (left panel) and 2D (right) modef^l. 
In the 3D model (left panel), non-axisymmetric structures are clearly seen. By performing 
a tracer-particle analysis, the maximum residency time of material in the gain region is 
shown to be longer for 3D than 2D due to the non-axisymmetric flow motions (see Figure [7]). 
This is one of advantageous aspects of 3D models to obtain the neutrino-driven explosions. 
On the other hand, our detailed analysis showed that convective matter motions below the 
gain radius beco me much more violent in 3D than in 2D, making the neutrino luminosity 
larger for 3D (see iTakiwaki et al.l (120121 ) for more details). Nevertheless the emitted neutrino 
energies are made smaller due to the enhanced cooling. Due to these competing ingredients, 
the neutrino-heating timescale becomes shorter for the 3D model, leading to a smaller net- 
heating rate compared to the corresponding 2D model (Figure El). Note here that the spectral 
IDSA scheme, by which the feedback between the mass accretion and the neutrino luminosity 
can be treated in a self-consistent manner (not like the light-bulb scheme assuming a constant 
luminosity), sounds quite efficient in the first-generation 3D simulations. 



As seen from Figure 0, an encouraging finding in ITakiwaki et al.l (120121 ) was that the 



15 The reason that the average shock of our 2D model expands much faster than theirs would come from the 
neglected effects in this work including general relativistic effects, inelastic neutrino-electron scattering, and 
cooling by heavy-lepton neutrinos. All of them could give a more optimistic condition to produce explosions. 
Apparently these ingredients should be appropriately implemented, which we hope to be practicable in the 
next-generation 3D simulations. 

16 Note that the polar axis is tilted (about n/4) both in the left and middle panel. 
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shock expansion tends to become more energetic for models with finer resolutions. These 
results would indicate whether these advantages for driving 3D explosions can or cannot 
overwhelm the disadvantages is sensitive to the employed numerical resolution^. To draw 
a robust conclusion, 3D simulations with much more higher numerical resolutions and also 
with more advanced treatment of neutrino transport as well as of gravity is needed, which 
could be hopefully practicable by utilizing forthcoming petaflops-class supercomputers. 



Fig. 6. — Volume rendering of entropy showing the blast morpholog y in our 3D (left) and 
2D (right) model for the 11.2M Q progenitor of I Woosley et al.l ( 120021 ) (at t = 178 ms after 
boun ce), respectively . The linear scale is indicated in each panel. These figures are taken 
from iTakiwaki et al.l (120 121 ). 



In addition to the 3D effects, impacts of GR on the neutrino-driven mechanism stand 
out among the biggest open questions in the supernova theory. It should be remembered that 
using newly derived Einstein equations (IMisner fc Sharplll964f ), the considerat i on of GR was 



standard i n the pioneering e ra of s upernova s imula tions (e.g., iMay &: White! ( 119661 )). One 
year after ICokate fc White! Jl966h . ISchwartzl Jl967h reported the first fully GR simulation 
of stellar collapse to study the supernova mechanism, who implemented a gray transport 
of neutrino diffusion in the ID GR hydrodynamicJ^I. Using GR Boltzmann equations de- 



17 It is of crucial i mportance to condu ct a convergence test in which a numerical gridding is changed in a 
systematic way (e.g. lHanke et al.l (|201ll )). 

18 Citing from his paper, "In this calculation, the neutrino luminosity of the core is found to be 10 54 erg/s, 
or 1/2 a solar rest mass per second !! .... This is the mechanism which the supernova explodes". The 




Fig. 7. — The left panel shows streamlines of selected tracer particles ad yecting through 



the s hock wave to the PNS seen from the polar direction in our 3D model (ITakiwaki et al. 



20121 ). Several surfaces of constant entropy marking the position of the shock wave (greenish 



outside) and the PNS (indicated by the central sphere) are shown with the linear scale given 
in the right bottom edge. The right panel shows the number of tracer-particles travelling in 
the gain region as a function of their individual residency time between the 2D and 3D model 
at 100 ms after bounce. If the left panel were for 2D models, the streamlines would be seen 
as a superposition of circles with different diameters. The maximum residency timescale for 
the 3D model i s shown to be longer t han that in the corresponding 2D model. These figures 
are taken from ITakiwaki et al.l (120121 ) . 
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Fig. 8. — Time evolution of neutrino-heating timescale (left) and total net rate of neutrino 
heating (right) i n our 2D (green line) and 3D model (red line). These figures are taken from 
Takiwaki et all (hoii ). 
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rived by iLindquistl ( 119661 ). IWilsonl (Il97l[ ) developed a ID GR-radiation-hydrodynamic code 



including a mor e realistic (at the time) description of the collisional term than the one in 



Schwartzl (119671 ). By performing ID GR hydrodynamic simulations that included a leakage 



van Riper & Lattimer 


1981; 


van Riper 


1982) 



pioneering studies, albeit using a much simplified neutrino physics than toda y, did provide 



a bott om-line of our current understanding of the supernova mechanism (see iBruenn et al. 
(200 l|) for a co mplete list of references for the early GR studies). In the middle of the 1980s, 
Bruennl (119851 ) developed a code that coupled ID GR hydrodynamics to the MGFLD trans- 
port of order (v/c) including the so-called standard set of neutrino interactions. Since the 
late 1990s, the ultimate ID simulations, in which the GR Boltzmann transport is coupled 
to ID GR hydrodynamic s , have been made feasible by Sumiyoshi-Yamada et al. (jYamada 

an d by Liebendorfer-Mezzacappa- 



19971 : lYamada etalllggd: ISumivoshi et al 



Bruenn et al. ( jMezzacappa fc Matzner 



2005b. 



1989; 



200 



Bruenn et al. 2001; Liebendorfer et al. 2001 



20041 ) (and by their collaborators). 



Among them, IBruenn et all ( 120011 ) firstly showed that the neutrino luminosity and the 
average neutrino energy of any neutrino flavor during the shock reheating phase increase 
when switching from Newtonian to GR hydrodynamics. They also pointed out that the 
increase is larger in magnitude compared to the decrease due to redshift effects and gravita- 



tional time dilation. By empl oying the current be st available weak interactions, iLentz et al. 



( 1201 ll ) reported the update of lBruenn et al.1 (j200lf ) very recently. They showed that the omis- 
sion of observer corrections in the transport equation particularly does harm to drive the 
neutrino-driven explosions. A disadvantageous trend of GR to produce the neutrino-driven 
explosions has been commonly observed in these full-fledged ID simulations; the residency 
time of material in the gain region becomes shorter due to the stronger gravitational pull. 
Due to this competition between the gain and loss effects in the end, GR works disadvanta- 
geously to facilitate the neutrino-driven explosions in ID. In fact, the maximum shock extent 
in the postbounce phase is shown to be 20% smaller when switching from Newtonian to GR 
hydrodynamics (e.g., Figure 2 in ILentz et al.1 (j201lf )). 



Among the most up-to-date multi-D models with spectral neutrino transport mentioned 
earlier (Table 1), the GR effects are at best attempted to be modeled by usin g a modified 



gravitational potential that takes into account a ID, post-Newtonian correction ( jBuras et al. 



neutrino luminosity rarely becomes so high in the modern simulations, but it is surprising that the potential 
impact of GR on the neutrino-heating mechanism was already indicated in the very first GR simulation. 

19 Very re cently, they reported their suc cess to develop the first multi-angle, multi-energy neutrino transport 
code in 3D ( Sumivoshi fc Yama da 2012). 
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2006 d ial: I Marek fc Jankall2009l ; iBruenn et al.ll20I0f ). A possible drawback of this prescription 
is that a conservation law for the total energy cannot be guaranteed by adding an artificial 
term in the Poisson equation. Since the energy reservoir of the supernova engines is the 
gravitational binding energy, any potential inaccuracies in the argument of gravity should 
be avoided. The re are a number of relativi stic sim ulations of massive stellar collapse in 
full GR (e.g., 2D Jshibata fe Sekiguchill2005al ) or 3D dshibata fc Sekiguchi|[2005bl lott et al. 
200 7dj), and references there i n) or using the conform a lly-fla tness approximation (CFC) (e.g., 
Dimmelmeier et al.l (J2002bJ); ICordero-Carrion et al.l (120091 ) ). Altho ugh extensive attem pts 
have been made to include micr ophysics such as by the Y e formula (iLiebendorferl l2005al ) or 
by the neutrino leakage scheme (jSekiguchill2010l ). the effects of neutrino heating have yet to 
be included in them, which has been a mai n hindrance to s tudy t he GR effects on the m ulti-D 
neutrino-driven mechanism (see, however, Mueller et al. ( 2012 ); Kuroda et al. ( 20121 )). 



Putting things together, a complete "realistic" supernova model should naturally be 
done in full GR(MHD) with multi-D GR Boltzmann neutrino transport, in which a mi- 
crophysical treatment of equation of state (EOS) and nuclear-neutrino interactions are im- 
plemented as realistically as possible. Unfortunately none of the currently published SN 
simulations satisfy the "ultimate" requirement (e.g., Table 1). In this sense, all the men- 
tioned studies employ some approximations (with different levels of sophistication) towards 
the final goal. 

In the same way, theoretical predictions of the SN multi-messengers that one can obtain 
by analyzing the currently available numerical results, cannot unambiguously give us the final 
answer yet. Again we hereby note the feature of this article which shows only a snapshot of 
the moving theoretical terrain. Keeping this caveat in mind, it is also true that a number 
of surprising GW features of CCSNe have been reported recently both by the first-principle 
simulations (e.g., in Table 1), and also by idealized simulations in which explosions are 
parametrically initiated mostly by the light-bulb scheme. As will be mentioned in the next 
section, the latter approach is also useful to get a better physical understanding of the GW 
signatures obtained in the first-principle simulations^!. Having shortly summarized a current 
status of CCSN simulations, we are now ready to move on to focus on the GW signatures 
from the next section. 



In this sense, these two approaches are complimentary in understanding the GW signatures. 
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2.2. Gravitational waves 



The paper by iMiillerl (119821 ) entitled as " Gravitational Radiation from Collapsing Ro- 



tating Stellar Cores" unquestionably opened our eyes to the importance o f making th e GW 



prediction based on realistic SN numerical modeling (see e.g., section 2 in lOttl (l2009b[ ) for a 
summary of more earlier work which had mainly focused on the GW emission in very ideal- 
ized systems such as in homogeneous spheroids and ellipsoids^. As one may expect from 
the title of his paper, rapid rotation, if it would exist in the precollapse iron core, leads to 
significant rotational flattening of the collapsing and bouncing core, which produces a time- 



Monchmever et al. ( 1991 ); Yamada fc Satol ( 1995); Izwerger &: Miiller ( 



(2003b 



( I2002a[ I2007L l2008h : IScheidegger et al 



2004bh: IShibata fc Sekiguchil (hoO^.lott et all feooi l2007iM: 





Miillei 


f 


1982f); 


(1997) 


; Kotake et al. 


; Dimmelmeier et al. 



( 120081 . l2010bf ) ) , and references therein) . 



As summarized by lOttl (l2009b[ ) , a number of important progresses have been recently 
made to underst and features of the bounce signals by ex tensive 2D GR studies using the CFC 
approximation (Dimmelmeier et al.ll2002al 120071 120081 ) and also by fully GR 3D simulations 
( jOtt et al.l 12004 ) both i ncluding realistic EO Ss and a deleptonization effect based on 1D- 
Boltzmann simulations ( iLiebendorferl l2005bl ) . 



For the bounce signals having a strong and characteristic signature, the iron core must 
rotate enough rapidly. However recent stellar evolution calculations suggest that rapid rota- 
tion assumed in most of the previous studies is not canonical for progenitors with neutron 
star formations dMaeder & Mevnetl I2OO0I ; Irleger et al] I2OO5I : btt et al.l l2006bh . To explain 



the observed rotation period s of radio pulsars, the precollapse rotation periods are estimated 
to be larger than ~ 100 sec ( jOtt et al.ll2006bl ). In such a slowly rotating case, the detection 



of the bounce signals bec omes very hard even by next-generation laser interferometers for a 
Galactic supernova (e.g., iKotake et al.l (j2004b[ )). 



Besides the rapid rotation, convective matter motions and anisotropic neutrino emission 
in the postbounce phase are expe cted to b e prim ary GW sources with comparable ampli- 



tudes to the bounce signals (e.g., iKotakd ( 1201 ll ) for a review). Thus far, various physical 



ingredients for producing asphericities and the resulting GWs in the post bounce phase have 



been studied, such as the roles of pre-c ollapse density inhomogeneities (IBurrows fc Haves 
1996 ; Miller &; Janka 1997 ; Fryer 2004bl laT). moderate rotation of the iron core ( Miiller et al. 



21 Needless to say, this kind of approach is still very important to extract the physics of the GW emission 
mechanism. 
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2004J), nonaxisymmetric rotational instabilities (iRampp et al.l 1 1998; Ott et al.l l2007bl ). g- 



mo des (lOtt et al.ll2006ai) and r-m o de s pulsations (lAndersson et al 



20101) of PNSs, and the 



SASI fcotake et all tooj l2009al lb}): iMarek et ali J2009h : iMurphv et al.l J2009[ )M Anion 



them, the most promising GW sources may be conv ection and SASI, because the degree of 
the initial inhomogeneities (lArnett fc Meakinll2011bl ) and the growth of rotational instabili- 
ties as well as the r-modes and g-modes pulsations are rather uncertain. 

Based on 2D simulations that parametrize the neutrino heating and cooling by the light- 
bulb scheme (which we shortly call as parametric SASI simulations hereafter), we pointed out 
that the GW amplitudes from anisotropic neutrino emission^] increase almost monotonically 
with ti me, and that such signals may be visible to next-generation detectors for a Galactic 
source (IKotake et al.ll2007l . l2009ai). By perf o rming such a parametric simulation but without 
the excision inside the PNS, IMurphv et al.l ( 120091 ) showed that the GW signals from matter 
motions can be a good indi cator of the exp l osion geometry. These features qualitatively agree 
with the ones obtained by lYakunin et al.l (120101 ) who reported exploding 2D simulations in 
which the ray- by-ray MGFLD neutri no transport is solv ed with the hydrodynamics (e.g., 
Oak-Ridge-I- simulations in Table 1). IMarek et al.l ( 120091 ) analyzed the GW emission based 
on their lo ng-term 2D ray-by-ray Boltzmann simulations, which seem very close to produce 
explosions (IMarek fc Jankal 120091 ) (e.g., MPA simulations in Table 1). They also confirmed 
that the GWs from neutr inos with continuou s ly grow i ng amplitudes ( but w ith the different 
sign of the amplitudes in IKotake et al.l (120071 . l2009al ) ; lYakunin et al.l ( 120 lOl ) ) , are dominant 
over the ones from matter motions. They proposed that the third-generation class detectors 
such as the Einstein Telescope are required for detecting the GW signals with a good signal- 
to-noise ratio. 



Regarding the GW predictions in 3D models, iMiiler fc Jankal (119971 ) coined the first 
study to analyze the GW signature of 3D non-radial matter motion and anisotropic neutrino 
emission from prompt convection in the outer layers of a PNS during the first 30 ms after 
bounce. Their first 3D calculations using the light-bulb recipe were forced to be performed 
in a wedge of opening angle of 60°. Albeit with this limitation (probably coming from the 
computer power at that time), they obtained important findings that because of smaller 
convective activities inside the cone with slower overturn velocities, the GW amplitudes of 
their 3D models are more than a factor of 10 smaller than those of the corresponding 2D 



2 see also lOtt et al.l (|2011l ) for the GW signals at the black hole formation. 



23 As anisotropic matter motions generate GWs, anisotrop i c neu t rino emission also g i ves ri se to GWs, 
which has been originally pointed out in late 1970's by Epstein ( 19781) : Turner fc Wagoner ( 1979) (see recent 
progress in Favata J 20101)'). I t is expected as a pr imary GW source a lso in gamma-ray bursts (|Hiramatsu et al 



20051 : ISuwa fc MuraselEoog ) and Pop III stars jSuwa et al.ll2007al ). 
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models, and the wave amplitudes from neutrinos are a factor of 10 larger than those due 
to non -radial matter mo t ions. With another pioneering (2D) study by Burrows & Hayes 
(1996) iBurrows fc Hayed (119961 ) , it is worth mentioning that those early studies had brought 
new blood into the conventional GW predictions, which illuminated the importance of the 
theoretical prediction of the neutrino GWs. 



A series of findings obtained by Fryer et al. in early 2000s (IFryer et al.l 12002k I2004J ) 
have illuminated also the importance of the 3D modeling. By running their 3D Newto- 
nian Smoothed-Particle-Hydrodynamic (SPH) code coupled to a gray flux-limited neutrino 
transport scheme, they studied the GW emission due to the inhomogeneous core-collapse, 
core rotation, low-modes convection, and anisotropic neutrino emission. Although the early 
shock-revival and the subsequent powerful explosions obtained in these SPH simulations 
have yet to be confirmed by other groups, their approach paying particular attention to the 
multiple interplay between the explosion dynamics, the GW signatures, the kick and spins 
of pulsars, and also the non-spherical explosive nucleosynthesis, blazed a new path on which 
current supernova studies are progressing. 

We also studied the GW signals from 3 D mod els that mimic neutrino-driven explosions 



aided by the SASI fcotake et al.l l2009al Jbl. | 201 ll ). In the series of our 3D experimental 



simulations, the light-bulb scheme was used to obtain explosions and the initial conditions 
were derived from a steady-state approximation of the postshock structure and the dynamics 
only outside an inner boundary at 50 km was solved. Based on the results, we show in the 
following that features of the gravitational waveforms obtained in 3D models are significantly 
different than those in 2D, which tells us a necessity of 3D modeling for a reliable prediction 
of the GW signals. 



2.2.1. Stochastic Nature of Gravitational Waves in 3D simulations 

Figure |H] shows the evolution of 3D hydrodynamic features from the onset of the non- 
linear regime of SASI (top left) until the shock break-out@(bottom right) with the gravita- 
tional waveform from neutrinos inserted in each panel. After about 100 ms, the deformation 
of the standing shock becomes remarkable marking the epoch when the SASI enters the non- 
linear regime (top left of Figure |9]). At the same time, the gravitational amplitudes begin 
to deviate from zero. Comparing the top two panels in Figure QUI which shows the total 
amplitudes (top panel, neutrino + matter) and the neutrino contribution only (bottom), it 



24 This corresponds to the shock emergence at the outer boundary of the computational domain (<~ 2000 
km in radius). 
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Fig. 9. — Four snapshots of the entropy distri butions of a r e presen tative 3D supernova 
explosion model (corresponding to model A in iKotake et al.l (l2009bl )). The second and 
fourth quadrant of each panel shows the surface of the standing shock wave. In the first and 
third quadrant, the profiles of the high entropy bubbles (colored by red) inside the section 
cut by the ZX plane are shown. The side length of each plot is 1000km. The insets show 
the gravitational waveforms from anisotropic neutrino emissions, with '+' on each curves 
representing the time of the snapshot. Note that the colors of the curves are taken to be the 
same as the top panel of Figure fTUl This figure is taken from IKotake et all (2009b). 
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Fig. 10. — Gravitational waveforms from neutrinos (bottom) and from the sum of neutrinos 
and matter motions (top), seen from the polar axis and along the equator (indicated by 
'Pole' and 'E quator') with p olariza tion (+ or x modes) for two representative 3D models of 
A and B (see iKotake et al.l (j2009bl ) for details). The distance to the SN is assumed to be 10 
kpc. 
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can be seen that the overall structures of the waveforms are predominantly determined by 
the neutrino-originated GWs with the slower temporal variations (> 30 — 50 ms), to which 
the GWs from matter motions with rapid temporal variations (< 10 ms) are superimposed. 

As seen from the top right through bottom left to right panels of Figure [9j the major 
axis of the growth of SASI is shown to be not aligned with the symmetric axis [:Z axis in 
the figure) and the flow inside the standing shock wave is not symmetric with respect to 
this major axis (see the first and third quadrant in Figure [H])- This is a generic feature in 
the computed 3D models, which is in contrast to the axisymmetric case. The GW ampli- 
tudes from SASI in 2D showed an increasin g trend with time due to the symmetry axis, 
along which SASI can develop preferentially f jKotake et al.l 120071 l2009af ). Free from such a 
restriction, a variety of the waveforms is shown to appear (see waveforms inserted in Figure 
5). Furthermore, the 3D standing shock can also oscillate in all directions, which leads to 
the smaller explosion anisotropy than 2D. With these two factors, the maximum amplitudes 
seen either from the equator or the pole becomes smaller than 2D. On the other hand, their 
sum in terms of the total radiated energy are found to be almost comparable between 2D 
and 3D models, which is likely to imply the energy equipartition with respect to the spatial 
dimensions. 

The (two) pair panels of Figure [10] show the gravitational waveforms for models with 
different neutrino luminosity. The input luminosity for the pair panels (models A (top two) 
and B (bottom two)) differs only 0.5^0. Despite the slight difference in the luminosities, 
the waveforms of each polarization are shown to exhibit no systematic similarity when seen 
from the pole or eq uator. This al s o refle cts a chaotic nature of the SASI influenced by small 
differences (see also iMiiller et al.l (j201lh ) . 



It should be noted that the approximations taken in the simulation, such as the excision 
inside the PNS with its fixed inner boundary and the light bulb approach with the isotropic 
luminosity constant with time, are the very first step to model the dynamics of the neutrino- 
heating explosion aided by SASI and study the resulting GWs. As already mentioned, the ex- 
cision of the central regi ons inside PNSs m ay hinder the efficient gravitational emis sion of the 



2007bl ; IScheidegger et a 



the PNSs ( Marek et al 



oscilla t ing neutron star (|Ott et al.ll2006af ) and the non-axisymmetric instabilities (lOtt et al. 



2010br) of the PNSs, and the enhanced neutrino emissions inside 
2009). Bearing these caveats in mind, a piece of encouraging news 



i s that the gravitatio nal waveforms obtained in the 2D radiation-hydrodynamic simulations 
( lYakunin et al.l 120101 ) are similar to the ones obtained in our 2D study using the light-bulb 



25 Although any seed perturbation would demonstrate the chaotic behavior, we simply chose the value of 
0.5 % as a reference. 
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scheme (IKotake et al.ll2009al ). More recently, iMuller et al.l (120111 ) c onfirmed the stochastic 



natur e by analyzing the GW features obtained in their 3D models (jWongwathanarat et al. 
20101 ). The obtained GW amplitudes as well as the degree of anisotropic neutrino emis- 



sion are almost comparable to our results, which are cons i derably smaller than 2D models 
rtMuller et allbooi iMarek et allbood : iMurphv et aDbood ; lYakunin et aDboioh . 
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Fig. 11. — Gravitational waveforms from the sum of neutrinos an d matter motions (le ft) and 
only from neutrinos (right) for a 3D model with rotation (from IKotake et al.l ( I201ll )). The 
time is measured from the epoch when the neutrino luminosity is injected from the surface 
of the neutrino sphere. For this 3D model with rotation, the rotational flow is imposed to 
advect to the PNS surface at around t = 400 ms. The supernova is assumed to be located 
at the distance of 10 kpc. 



2.2.2. Breaking of the Stochasticity Due to Stellar Rotation 



More recently, we st udied the effects o: 



GWs mentioned above (IKotake et al.l 12011 



stellar rotation on the stochastic nature of the 
). In 3D, t he m odes of SASI are divided into 
( 120081 ) ). Asymmetric m = modes 



Iwakami et al 



sloshing modes and spiral modes (e.g., 
so far studied in 2D models and axisymmetric m/fl modes are classified into the sloshing 
modes, where m stands for the azimuthal index of the spherical harmonics Y™. In the latter 
situation, the ±m modes degenerate so that the +m modes has the same amplitudes as 
the — m modes. If random perturbations or uniformly rotating flows are imposed on these 
axisymmetric flows in the postbounce phase, the degeneracy is broken and the rotational 
modes emerge (IBlondin fc Mezzacappall2007l ; llwakami et al.ll2009l ). In this situation, the +m 
modes has the different amplitudes from — m modes. Such rotating non-axisymmetric m/fl 
modes are called as spiral modes. These non-axisymmetric modes are expected to bring 
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Fig. 12. — Partial cutaway of the entropy isosurfaces and the velocity vectors on the cutting 
plane for a 3D model that includes rotation. Left and right panels are for the equatorial and 
polar observer, respectively. The insets show the gravitational waveforms with '+' on each 
curves representing the time of the snapshot. Note that the colors of the curves are taken to 
be the same as the top panel of Figure [TUJ This figure is taken from lKotake et al.l ( 120111 ). 
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about a breakthrough in our supernova theory, because they can help to produce explo sions 
more easily compared to 2D due to its extra degree of freedom (INordhaus et al.l 12010) , and 



also b e cause they may have a potenti al to generate pulsar spins (IBlondi n &: Mezzacappa 



f j2007h : iRantsiou et al l fcoiol ) , see also lYamasaki fc Foglizzol (j2008l ) ; IWongwathanarat et al.l 



( 201ol )). These finding illuminating the importance of stellar rotation motivated us to clarify 
how rotation that gives a special direction to the system (i.e., the spin axis), could affect the 
stochastic GW features which we observed in the absence of rotation (section |2.2. ip . 

Figure [TT] shows the gravitational waveform for a typical 3D model with rotation (left:total 
amplitudes, right: neutrino only). To construct a model with rotation, we giv e a uniform ro- 



tation on the flow advecting from the outer boundary of the iron core as in llwakami et al. 



( 120091 ). whose specific angular momentum is assum ed to agree with recent stellar evolution 
models (IMaeder &: Meynetl I200CH : iHeger et al.l 120051 ) . Comparing to Figure [TU] (for models 
without rotation), one can clearly see a sudden rise in the GW amplitude after around 500 
ms for the rotating model (blue line in Figure [TT]) . which is the plus mode of the neutrino GW 
seen from the equator. By systematically changing the initial angular mome ntum and the 



i nput neutrino luminosities from the PNS, we computed fifteen 3D models in iKotake et al. 



( ]201ll ) and found that the GW features mentioned above were common among the models. 



Figure [12] illustrates a typical snapshot of the flow fields for the rotating model (cor- 
responding to the one in Figure [TT]) when the spiral SASI modes have already entered the 
non- linear regime, seen from the pole (right panel) or from the equator (left panel), respec- 
tively. From the left panel, one may guess the presence of the sloshing modes that happen 
to develop along the rotational axis (z-axis) at this epoch. It should be emphasized that 
although the dominance of /i^ q + observed in the curre nt 3D simulations is similar to the one 
obtained in previous 2D studies (IKotake et al.ll2009al ). its origin has nothing to do with the 
coordinate symmetry axis. The preferred direction here is determined by the spin axis. Free 
from the 2D axis effects, the major axis of the sloshing SASI mode changes stochastically 
with time, and the flow patters behind the standing shock simultaneously change in every 
direction like the non-rotating models. As a result, the sloshing modes can make only a 
small contribution to the GW emission. The remaining possibility is that the spiral flows 
seen in the right panel should be a key importance to understand the GW feature mentioned 
above. In fact, by analyzing the matter distribution on the equatorial plane, we find that 
the compression of matter is more enhanced in the vicinity of the equatorial plane due to 
the growth of the spiral SASI modes, leading to the formation of the spiral flows circulat- 
ing around the spin axis with higher temperatures. As a result, the neutrino emission seen 
parallel to the spin axis becomes higher than the ones seen from the other direction. Remem- 
berin g that the lateral-ang le (9) dependent function of the GW formulae (e.g., in equation 
(9) in lKotake et al.l (l2009al )) is positive near the north and south polar caps, the dominance 
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Fig. 13. — Spectral distributions of GWs from matter motions ("Matter") and neu- 
trino emission ("Neutrino") s een from th e pole or the equator for a representa- 
tive 3D rotating model (e.g.. iKotake et al.l (120111)') with the expected detection lim- 
its of 
LIGO 



AMA30 (lAndo & the TAMA Collaboration! 120051 ). first LIGO and advanced 



Thornd Il995l). Large-scale Cryogenic Gravitational wave T elescope (LCGT ) 



flKuroda fc LCGT Collaboration! BOlOh and Fabry-Perot type DECIGO (iKawamural 120061 ). 
It is noted that h r h»r is the characteristic gravitational wave strain defined in 



Flanagan fc Hughes! (119981 ). The distance to the supernova is assumed to be 10 kpc. Note 
that for the matte r signal, the + mode seen from the polar direction is plotted (from 
Kotake et all (boilh ). 
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of the polar neutrino luminosities leads to make the positively growing feature of h e J^ in 
Figure [TT] (blue line). From the spectral analysis of the gravitational waveform (Figure [T2l . 
it can be readily seen that it is not easy to detect these neutrino-originated GW signatures 
with slower temporal evolution (> O(lO)ms) by ground-based detectors whose sensitivity is 
limited mainly by the seismic noises at such lower frequencies. However these signals may 
be detecta ble by the recentl y proposed future space interferometers like Fabry-Perot type 
DECIGO (jKawamural (120061 ). black line in Figure [T2|) . Contributed by the neutrino GWs 
in the lower frequency domains, the total GW spectrum tends to become rather flat over 
a broad frequency range below ~ 100 Hz. These GW features obtained in the context of 
the SASI-aided neutrino-driven mechanism are different fro m the ones expected i n the other 



candidate mechanisms, s uch as the MHD mechanism (e.g.. iTakiwaki fc Kotakd (120111 )) and 



the acoustic mechanism (jOtt et al.ll2006al ). Therefore the detection of such signals could be 



expected to provide an important probe into the explosion mechanism (e.g., Ott (j2009bl Jah ) . 



We like to draw a caution that most of the 3D models cut out the PNS and the neutrino 



transport is appro ximated by a simpl e light-bulb scheme (IKotake et al.ll201ll ) or by the gray 
transport scheme (jMuller et al.ll201ll ). Needless to say, these exploratory approaches are but 
the very first step to model the neutrino-heating explosion and to study the resulting GWs. 
As already mentioned, the excision of the central regions inside PNSs truncates the feedback 
between the mass accretion to the PNS and the resulting neutrino luminosity, which should 
affect the features o f the neutrino GW s. By the cut-out, efficient GW em ission of the oscil- 
lating neutr on star (jOtt et al.ll2006al ) and non-axisymmetric instabilities (IScheidegger et al. 



2008 



s 



201QbJ) of the PNSs, and the enhanced neutrino emissions inside the PNSs (IMarek et al.l 



20091 ) cannot be treated in principle. To elucidate the GW signatures in a more quantitative 
manner, full 3D simulations with spectral neutrino transport are apparently needed (e.g., 
section 2.1). This is unquestionably a vast virgin territory awaited to be explored for the 
future. 



2.3. Explosive nucleosynthesis 

In this section, we proceed to discuss possible signatures of supernova nucleosynthe- 
sis. The study of nucleosynthesis is of primary importance to unveil the origins of heavy 
elements. It could also provide a valuable information of the ejecta morphology by observ- 
ing the aspherical distributions of the synthesized elements especially for a nearby CCSN 



-34- 



eventl I. In the following, we first present a short overview paying particular attention to 
explosive nucleosynthesis, and then discuss possible observational signatures that would im- 
print information of multidimensionalities of the supernova engine. 

When in a successful explosion the shock passes through the out er shells, its high tem - 



perature induces an explo si ve nucleosynthesis on short timescales (e.g.. lWooslev et a 



(2002); 



Woosley fc Weaverl (119951 ) ; iHashimotd (119951 ). and collective references in lJerkstrandl ( 120111 )). 
The observational determination of the masses of the three main radioactive isotopes 56 Ni, 
57 Ni, and 44 Ti sets one of the main constraints on the explosion dynamics, because the 
production of these elements is sensitive to the track of d ensity and temperature that the 
expanding material traces (e.g.. IWoosley Hoffman! (Il99lf )). During the shock propagation, 



iron group elements such as 56 Ni and its daughter nucleus 56 Co are predominantly produced, 
which are radioactive with a life-time of 8.8 days and 111.5 days, respectively. Most CCSNe 
enter the so-called nebular phase after the first few months when the expanding ejecta be- 
comes optically thin in the continuum. In the early nebular phase, 56 Co is the major nuclear 
power source. As long as the decay particles are trapped by the ejecta, the radiation energy 
supplied by radioactivity is emitted instantaneously, so that the light curve can be described 
by an exponential decay with time, simply tracing the decay of the 56 Co nuclide. To explain 
the bolo metric light curve of SN1987A in such a phase, the 56 Ni mass was determined to be 
0.07M© faouchet et alill99lh . 



After several years of explosion, the radioactive output from the ejecta no longer bal- 
ances with the inst antaneous input by radio activity, because the reprocessing timescale is 
going to be longer (IFransson fc Kozmalll993l ). The bolometric light curve is affected by the 
delayed release of the ionization energy. After that, a self-consistent modeling is needed, 
in which one should include a detailed calculation of the gamma-ray /positron thermaliza- 
tion and a determination of the time-dependent t emperature, ionizati o n, and excitation ( e.g., 



Fransson fc Kozmal (119931 1 ; lMazzali et all ( 12001f ); lMaurer et al I (I201lh: 



Jerkstrand et al 



(1201 if ) 



and references therein). Such a time-dependent modeling by lFransson &: Kozmal (119931 ) re- 
vealed the 57 Ni mas s of~ 3.3 x 10~ 3 M(p of SN1987A, which agrees well with observations (e.g., 
Varani et al. ( 1990 ); Kurfess et al. ( 1992 ) and collective references in Filippenko (jl997 )). 



By the similar reason to 57 Ni just mentioned above, t he de termination of the 44 Ti is also 
complicated. The most recent study by I Jerkstrand et al.l (120111 ) gives an estimate of the 44 Ti 
to be 1.5^q;5 x 10 ~ 4 M(^, which i s in go od agree ment with the eight-yea r spectrum analysis 
of SN1987A (e.g., IChugai et al.l (Il997| ). see also iLundqvist et al.l ( 120011 )). As shown above, 



26 Note that nucleosynthesis is not critical for the modeling of the light-curve and spectra for the most 
frequent types of SNe II-P. 
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the amount of 44 Ti is typically one order-of-magnitude smaller than that of 57 Ni, however it 
is crucially important for young supernova remnants due to its long life-time (~ 86 years). 
It is worth mentioning that NASA will launch the satellite NuSTAR (Nuclear Spectroscopic 
Telescope Array) to study 44 Ti production in CCSNe. The detector will be able to map out 
the 44 Ti distribution of the supernova remnant Cassiopeia A and can get velocity distribu- 
tions of the 44 Ti in SN 1987A. By comp aring detailed mod e ling of the SN nuc l eosynt hesis in 
the context of 2D and 3D models (e.g., iFryer et al.l (120061 ); iMagkotsios et al.l ( 120101 )). these 
are expected to both provide direct probes of the explosion asymmetry. 



Ever since SN1987A, challenges to the class i cal spherical modelin g (jHashimotol 11995 
Woosley fc Weaver! Il995l ; iThielemann et al.lll996t iRauscher et al.l 120021 ) have been building 
also in the SN nucleosynthesis (likewise in the explosion theory and GWs mentioned so far). 
For many years it has been customary to simulate explosions and the effects of the shock wave 
on the explosive nucleosynthesis by igniting a thermal bomb in the star's interior or by initi- 
ating the explosion by a strong push with a piston. 2D simulations with manually imparted 
asymmetries showed that bipolar explosion scenar ios could account for e nhanced 44 Ti synthe- 



sis along the poles as indicated in SN1987A ( e.g.. iNagataki et al.l (119971) '). More recently, 3D 



effects have been more elaborately studied ( jHungerford et al.l (120051 ) . see also lYoung et al. 



(120061)) as well as the im pacts of different explosions by em ploying a number of progeni- 



tors (IJoggerst et al.l 120101 ) or by assuming a jet-like explosion (ICouch et al 



20091 ). which is one of the possible candidates of hypernovae 



(Couch et al. 


2009; 


Tominaga 




Maeda & Nomoto 


(2003 





In addition to the above-mentioned work, nucleosynthesis in a more realistic simula- 
tion that mode l s mult i -D ne u trino-driven explosi o ns ha s been also extensively studied (e.g., 
Kifonidis et al.l ( 120031 . 120061 ) ; iGawryszczak et al.l (120101 ) and references therein). Although 
a small network has ever been included in the computations, these 2D simulations employ- 
i ng a light-bulb scheme dKifonidis et al. 2003 ) or a more accurate gray transport scheme 
( iScheck et al.ll2006l ; iKifonidis et al.l 120061 ) have made it possible to elucidate nucleosynthesis 
inside from the i ron co re after the shock-revival up to explosion in a more consistent manner. 
Kifonidis et al.l ( 120061 ) demonstrated that the SASI-aided low-mode explosions can naturally 
explain masses an d distributions of the synthesized elements observed in SN1987A. Their re- 
cent 3D results by lHammer et al.l ( 120101 ) show that the 3D effects change the velocity profiles 
as well as the growth of the Rayleigh- Taylor instability, which affects properties of the SN 
ejecta. In simulations with spectral neutrino transport, a new nucleosynthesis process, the 
so-called z/p process, is reported to successfu l ly exp l ain some light proton - rich(p-) nuclei in- 
cludin g 92 < 94 Mo and 96 ' 98 Ru (e.g JPruet et al.l J20061 ) ; iFrohlich et all fco06h ; lwanaio & Janka 
(120 111 ) and references therein). It should be noted that self-consistent simulations are cur- 
rently too computationally expensive to follow the dynamics of the expanding shock far 
outside the central iron core 1 s after bounce) where explosive nucleosynthesis takes 
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place. On the other hand, the light-bulb scheme is not accurate to determine a sensitive 
balance between neutrino captures proceeding via v e or u e , which is decisive for quantifying 
the processes. Therefore the two approaches, namely light-bulb scheme vs. self-consistent 
neutrino transport, may be regarded as playing a complimentary role at present. 

In the next section, we briefly summarize our findings on explosive nucleosynthesis in our 



2D models that utilize the light-bulb scheme to trigger explosions (e.g.. lFujimoto et al.l (120111 ) 
for more details). By changing neutrino luminosities from PNSs systematically (L v ), we 
discuss how the multidimensionality formed by the SASI and convection could impact on the 
explosive nucleosynthesis. We employ a non-rotating 15M Q star with solar metalicity, which 
has been a best-studied supernova progenitor. The abundance pattern of the synthesized 
elements is estimated by a post-processing procedure, in which we have adopted a large 
nuclear network continuousl y maintained by th e Joint Institute for Nuclear Astrophysics 



(JIN A) REACLIB project (jCyburt et al.l l2010l ). Our readers might wonder why we are 



returning to 2D results here again after referring to 3D studies in section 12.21 This is simply 
because our 3D project for the nucleosynthesis is currently in progress, and we like to note 
again the feature of this review article which shows only a snapshot of the moving theoretical 
terrain. 



2.3.1. Symmetry breaking in explosive nucleosynthesis 

Figure [TH shows how explosion energies (top panel), explosion timescales as well as 
the PNS masses (bottom panel) chang e with the inpu t neut rino luminosity (L Ve ) for the 



2D parametric explosion models (e.g., iFujimoto et al.l ( 1201 f I ) for more details). As seen 



higher L Ue makes explosion energy greater, the onset of explosion (t exp ) earlier, the PNS 
masses (Mpns) smaller. The PNS masses range in (1.54 — 1.62)M for models with higher 
luminosity (L Ue > 4.25 x 10 52 erg s -1 , panels (b), blue line), which a re comparable to th e 



baryonic mass (around 1.5M ) of neutron stars observed in binaries (ISchwab et al.ll2010l ). 
The explosion energies for these high luminosity models (panel (a), red line) are also close 
to the canonical supernova kinetic energy of 10 51 erg. To discuss explosive nucleosynthesis 
in the following, we thus choose to focus on the model with L Ve = 4.5 x 10 52 ergs _1 as a 
reference. 

For the model, Figure [TBI shows a snapshot of internal energy (left panel) and abundances 
of ejecta (O, Si, Fe) when the SASI-aided low-mode shock is propagating in the O-rich layers 
at around ~ 20000 km along the pole. Ejecta with higher internal energies (left panel, colored 
by red) concentrate on the shock front particularly in polar regions (6 < n/4,6 > 37r/4), 
where both Si and O burning proceed to produce abundantly 56 Ni and 28 Si. These aspherical 
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Fig. 14. — Panel (a) shows explosion energies (red line) vs. input neutrino luminosities {L Ve ), 
where green and blue line shows contributions from kinetic and thermal energy, respectively. 
Panel (b) shows the mass of protoneutron star (M PNS ) and the explosion timescales (t cxp ) 
as a function of L Ue . M PNS is estimated as the enclosed mass exceeding 10 12 g cm -3 and t exp 
as the time scale when the mass ejectio n rate at 100km drops down to O.lM /s in our 2D 
simulations (from lFujimoto et al.l (120111 ) ) . 
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element distributions mi ght be respons i ble fo r anisotropies o f the SN ejecta as observed in 
SN1987A. More recently, IrCimura et al.l (120091 ): llJchida et al.l (120091 ) have analyzed the metal 
distribution of the Cygnus loop and pointed out that the progenitor of the Cygnus loop is a 
CCSN explosion whose progenitor mass ranges in ~ 12 — 15M Q . Since the material in this 
middle-aged supernova remnant has not been completely mixed yet, the observed asymmetry 
is considered to still remain a trace of inhomogeneity produced at the moment of explosion. 
The asymmetries of the SN ejecta obtained in our 2D explosi on model of a 15 Mq p rogenitor 
have a close correlation with the ones in the Cygnus loop (see iFujimoto et al.l ( 120111 ) for more 
detailed comparison). This might allow us to speculate that globally asymmetric explosions 
induced by the SASI-aided neutrino-heating mechanism could account for the origin. 

Concerning the synthesized amount of 56 Ni and 44 Ti, they are ~ O.O6M and 4x 1O _5 M , 
respectively for the model (left panel in Figure [ToT). For SN1987A, m asses of 56 Ni and 44 Ti 
are deduced to be ~ O.O7M Jshigevama et al.l Il988l : IWooslevlll988h and 1 - 2 x 1O~ 4 M 



(lNagatakill2000l . and references therein). The obtained mass of 44 Ti is, theref ore, not enough 
for SN1987A as well as for Cas A (1.6±gJ x 1O" 4 M ) jRenaud et a lj|2006h. The shortage 
of 44 Ti is a long-standing problem, which might be solved by 3D effects (IHammer et al. 



20101 ). For the abundance patterns, our reference model (top panel in Figure [TBI) is shown to 
reproduce a similar trend to the solar system (right panel, top), which is in sharp contrast 
to the model with lower neutrino luminosity (right panel, bottom). 

Finally we note that there should are at least two barriers to get over in the current- 
generation studies of the CCSN nucleosynthesis. As mentioned already, the first one is the 
need of radiation hydrodynamic simulations that follow the dynamics from gravitational 
collapse, shock-revival, shock propagation in the stellar mantle, to stellar explosions in a 
self-consistent manner. The second one is the need of computing light curves, spectra, 
and the degree of polarizations in 2D or 3D hydrodynamic models, for which one needs to 
solve the multi-D photon transport coupled with multi-D hydrodynamics. Unfortunately 
however, the marriage of the two items, albeit ultimately needed for clarifying the photon 
messengers, may not be done immediately due to their computational expensiveness. For 
example, one needs to follow the dynamics more than ~ 1 day after the onset of gravitational 
core-collapse for computing the explosive nucleosynthesis, because the envelop of a typical 
supernova progenitor extends up to a radius of ~ 10 12 cm. It takes furthermore more than 
~ 1 week before the shock propagation enters to the so-called homologou s phase. At present. 



the be st available numerical simulation to this end is limited to 2D (e.g.. iGawryszczak et al. 



(120101 )). For discussing anisotropies in emission-line profiles, one needs to follow the dynamics 
later than ~ 1 year after explosions when the remnant becomes transparent. It is a very 
challenging (but very important) task for the first principle simulations to overcome the 
big gaps regarding the very different timescales. A encouraging news is that several groups 
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Fig. 15. — Distribution of internal energy (left panel) and abundances of the ejecta (0, Si, 
Fe from left to right, normalized to solar) when the non-spherical shock is propagating in the 
O-rich layers at around ~ 20000 km along the pole (~ 1.5 s after bounce). For this model, 
the input neutrino luminosity is set as L Ve = 4.5 x lO^ergs" 1 . N ote that Fe in the right 
panel is the sum of 56 Fe, 56 Ni, and 54 Fe (from iFujimoto et al.l (l201lh ) . 
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Fig. 16. — Left panel shows masses of ejecta (M^in), and masses of radioactives ( 56 Ni and 
44 Ti). Right panel shows element abundance yields (Xi/X Q; i, normalized to solar) vs. mass 
number for our reference model of L Ue = 4.5 x 10 52 erg s _1 and for (b) L Ve = 3.9 x 10 52 erg s~ x 
Thick horizontal-dashed lines represent a factor equals to that of 16 0, while two normal and 
two thin lines denote a fact or equals to that of 16 Q times 2, 1/2, 3, and 1/3, respectively. 
These figures are taken from lFujimoto et al.l (120111 ) . 
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are pursuing it with the use of advanced n umerical techniques suc h as the adaptive-mesh- 



refine ment approach (INordhaus et al.ll2010l ) or Ying-Yang gridding (jWongwathanarat et al. 



2010l ) . With an accelerating power (like peta- or exa-scale) of supercomputers, it would be 



unsurprising that the next (or the next after next!) generation supernova modellers can 
execute the ultimate ab-initio simulations, whose findings would be immediately used to 
make a more precise prediction of the photon messengers. 



3. MHD Mechanism 



In this section, we move on to focus on the MHD mechanism. After we briefly summarize 
the current status of this topic below, we mention the explosion dynamics that proceeds 
by the field-wrapping mechanism in section 13.0.21 and discuss properties of the so-called 
magnetorotational instability (MRI) in sections 13.0.31 and 13.0.41 based on our preliminary 
results. Possible signatures of GWs and neutrinos are discussed in section 13.11 and 13.21 
respectively. 

Numerical simulations of MHD stellar explosions have started already in early 1970's 



1976 



shortly after the discovery of pulsars (jLeBlanc &; Wilsonll 19701 ; lBisnovatyi-Koganlll970t iBisnovatyi-Kogan et 



Miiller fc Hillebrandtlll979l ; ISymbalistylll984j ). However, it is rather only recently that 



the MHD studies come back to the front-end to pics in the supernova research followed 



by a number of extensive A 


m 


D simulations ( 


e.g., Ardelif 


in et a 


. (2000 


); Yamada & Sawa 


i 


(2004) 


: Kotake et al. 


]2004t 


a, 


2005); 


Obergau 


ineer et al. 


(2006a 


b); 


Burrows et al. 


(2007a) 




Cerda-Duran et al. (I 


2007); i 


Miwa et al 


( 


2007b « 


3); Komissarov & '. 


Sarkovl ( 


2007):lTakiwaki et al. 


(2004, 


2009): Scheideeeer et al. 


(2010a 


); 


Obereaulineer & Janka ( 


2011 


) for references therein) . 



Main reaso ns for this activity are observations in dicating very asymmetric exp 



2001 



2002 ). and the interpretation of magnetars (jDuncan Sz Thompson 1992 ; Lattimer fc Prakash] 



2007 ). collapsars and their reley ance to gamma-ray bursts ( Dessart et al. ( 2008^ ; Woosley fc He 



osions (IWang et al. 



d2006h : lYoon fc Langej(|2005bh ) as a possible outcome of the magnetorotational core-collapse 
of massive stars. 

The MHD mechanism of stellar explosions relies on the extraction of rotational free 
energy of collapsing progenitor core via magnetic fields. Hen ce a high angular momentum 
of the core is preconditioned for facilitating the mechanism ( jMeier et al.lll976l ). Given (a 
rapid) rotation of the precollapse core, there are at least two ways to amplify the initial 
magnetic fields to a dynamically important strength, namely by the field wrapping by means 
of differential rotation that n aturally develops in the collapsing core, and by the MRI (MRI, 



^gcr 



sec 



Balbus fe Hawleyi (Il998[ )). Each of them we are going to review briefly in the following. 



Fig. 17. — Three dimensional plots of entropy with the magnetic field lines (left) and the 
streamlines of the matter (right) during the j et propagation a t 20 m s and 94 ms after bounce, 
respectively (for model B12TW1 taken from Takiwaki et al. ( 20091 )). The outer edge of the 
sphere colored by blue represents the radius of 7.5 x 10 7 cm. These panels highlight not only 
the wound up magnetic field around the rotational axis (left), but also the fallback of the 
matter from the jet head of the downwards to the equator, making a cocoon-like structure 
behind the jet (right). 
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Fig. 18. — Left and middle panels show magnetic pressure (red line) vs. ram pressure (blue 
line) for a typical MHD model (same as Figure ITT1) along the polar axis (left panel) or the 
equatorial plane (middle panel) at 2 ms after the revival of the stalled bounce shock. The 
matter pressure is shown by green line as a reference. The right panel shows a velocity profile 
along the pole measured from the shock-stall (0.0 ms in the panel). For the equator, the 
magnetic pressure is much less than the ram pressure (middle panel), while the magnetic 
pressure amplified by the field wrapping along the pole becomes as high as the ram pressure 
of the infalling material at the shock front (left panel, note that the shock position can be 
inferred by the discontinuity at around 150 km), leading to the MHD-d riven shock formation 
(right panel). These figures are taken from iTakiwaki &: Kotakd (120111 ). 



3.0.2. Field-wrapping mechanism 
Three dimensional plots of Figure IT71 are useful to see how the field wrap ping occurs. For 



the 2 D model taken from our 2D special relativistic (SR) MHD simulation (ITakiwaki et al. 



20091 ). the precollapse magnetic field is set to be as high as 10 12 G that is uniform and parallel 
to the rotational axis and the initial (5 parameter (ratio of initial rotational energy to the 
absolute value of the initial gravitational energy) is taken to be 0.1 % with a uniform rotation 
imposed in the iron core. From the left panel, it can be seen that the field lines are strongly 
twisted around the rotational axis. From the induction equation of i deal MHD equation s, 



the time evolution of the toroidal fields (B^) can be expressed as (e.g.. iMeier et al.l ( 119761 )) 



at ~ B *K X 6xh (1) 

where X denotes the distance from the rotational axis and Bx represents the X component 
of the poloidal fields in the cylindrical coordinates (X, Z). Given that a precollapse iron core 
has strong magnetic field (such as Bx,o ~ 10 12 G with Bx,o representing the initial poloidal 
fields) and rapid rotation (such as P^t ^ 4 s with P ini t being the precollapse rotation period), 
the typical amplification of B^ near core bounce may be estimated as 

AB ^ M >* G (^)(^)(Jiy, (2) 
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where t represents the typical timescale when the field amplification via the wrapping pro- 
cesses saturates, P TOt denotes the rotational period near bounce. Equation (2) shows that 
the toroidal fields grow lin early with time by wrappi ng the poloidal fields (Bx), which is the 



so-called Q dynamo (e.g.. ICerda-Duran et al.l (120071 )). The 10 15 G class magnetic fields are 



already dynamically important strength to affect the postbounce hydrodynamics as will be 
mentioned below. 

The white lines in the right panel in Figure [T71 show the streamlines of matter. A fallback 
of material just outside of the jet head advecting downwards to the equator (like a cocoon) 
is seen. In this jet with a cocoon-like structure, the magnetic pressure is always dominant 
over the matter pressure. As estimated above, the typical field strength behind the shock is 
~ 10 15 G. The critical strength of the toroidal magnetic field to induce the magnetic shock 
revival (i.e. the revival of the stalled bounce shock due to the MHD mechanism) is estimated 
as follows. The matter in the regions behind the stalled shock is pushed inwards by the ram 
pressure of the accreting matter. This ram pressure is estimated as, 



P = 4 x 10 



2.s 



P 



Av 



erg/ cm 3 



(3) 



10 10 g/ cm3 / V 2 x 10 9 cm/s 
where p and Av denotes typical density and radial velocity in the vicinity of the stalled shock. 
When the toroidal magnetic fields are amplified as large as ~ 10 15 G due to the field wrap- 
ping behind the shock, the resulting magnetic pressure, |^ ~ (10 15 ) 2 /87r ~ 10 29 erg/cm 3 , 
can overwhelm the ram pressure, leading to the magnetic shock revival (e.g., Figure [T8|) . 
The onset timescale of the magnet ic shock reviva l is shar ply dependent on the precollapse 
magnetic fields and rotation (e.g., iBurrows et al.l (j2007al )). As the initial field strength is 
larger with more rapid rotation imposed, the interval from the shock stall to the MHD-driven 
shock revival becomes shorter. The speed of the jet head is typically mildly relativistic (at 
most ~ 0.4c), with c being the speed of light (see also the right panel of Figure [T8|) . At the 
shock breakout of the iron core, the explosion energies generally exceed 10 51 erg for models 
that produce MHD explosions in a shorter delay after the stall of the bounce shock. These 
features are in accord with those obtained in more detailed MHD simulations with spectral 
neutrino transport (e.g., IBurrows et al.l (j2007al ); iDessart et al. 



3.0.3. On the MRI 



We are now in a position to discuss possible impacts of the MRI in supernova cores. 
Akiyama et al.l ( 120031 ) was the first to point out that the interfaces surrounding the nascent 



27 See, e.g., Takiwaki fc Kotakel ( 2011 ) for a more detailed comparison. 
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PNSs quite generally satisfy the MRI instability criteria. Therefore any seed magnetic fields 
can be amplified exponentially in the differentially rotating layers, much faster than the 
linear amplification due to the field wrapping (e.g., Equation After the MRI enters to 
the saturated state, the field strength might reach ~ 10 15 ~ 16 G, which is high enough to affect 
the supernova dynamics. Not only in the field amplificat i on, th e MRI plays a cruc i al rol e 
also in operat i ng th e MHD turbulence (see lHawley et al.l (119951 ); iBalbus fc Hawleyl (119981 ); 
Masada et al.l (120061 )). The turbulent viscosity sustained by the M RI can convert a f r action 
of the shear rotational energy to the thermal energy of the system. iThompson et al.l (120051 ) 
suggested that the additional energy input by the viscous heating can help the neutrino- 
driven supernova explosion. Followed by the exponential field amplification and additional 
heating, a natural outcome of the magnetorotational core-collapse is expected to be the 
formation of energetic bipolar explosions, which might be observed as hypernovae (e.g., 
section |2~51) . 

Note again that bipolar explosions obtained in the previous section with the assumption 
of a strong precollapse magnetic field (> 10 12 G), a r e pre dominantly driven by the field 
wrapping processes, not by the MRI. IShibata et al.l (120061 ) pointed out in their fully 2D 
GRMHD core-collapse simulations that more t han 10 grids a r e at l east needed to capture 
the fastest-growing mode of the MRI (see also lEtienne et al.l (120061 )). As well-known, the 
growth rate of MRI-unstabl e modes depends on the product of the field strength and the 
wave number of the mode (IBalbus &; Hawleyl Il998l ). When a precollapse rapidly core is 
strongly magnetized as 10 11 G (like in the case of collapsar progenitor), the wavelength is 
on the order of km in the postbounce core. In this case, the ex ponential field grow th of 
the MRI was successfully captured by high-resolution simulations ( Shibata et ali 2006). On 



the o ther hand, the fastest growing modes drop below several meters (jObergaulinger et al. 



20091 ) for canonical supernova progenitors, which rotates much more s lowly with weaker 
initial fields (~ 10 9 G) as predicted by recent stellar evolution models (IMaeder fc Meynet 



2000l ; iHeger et al.ll2005l ). At present, it is computationally too expensive to resolve those 



small scales in the global MHD simulations, typically more than two or three orders-of- 
magnitudes smaller than their typical finest grid size. To reveal the nature of the MRI, 
local simulations focusing on a small part of the MRI-unstable regions a re expected to be 
quite useful as traditionally studied in the context of accretion disks (see IBalbus Hawley 
jl998h ). 



Obergaulinger et al.l (120091 ) were the first to study the growth and saturation level of 
the MRI in the supernova environment^. To ease a drawback of the local shearing box 



2.s 



sec, 



Zhang et al.l (j2009h for extensive simulations that focus on the MRI in relativistic outflows in the 



context of gamma-ray bursts and active galactic nuclei. 
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simulation, they employed the shearing disk boundary conditions by which global radial 
density stratification can be taken into account. By performing such a semi-global simulation 
systematical ly in 2D and 3D , they derived scaling laws for the termination of the MRI. As 
estimated in lAkiyama et al.l ( 120031 ). the MRI was shown to amplify the seed fields exceeding 
10 15 G. These important findings may open several questions that motivate us to join in 
this effort, such as how the nonlinear properties as well as the scaling laws could be in the 
subsequent non-linear state and whether the viscous and resistive heating driven by the MRI 
turbulence could or could not affect the supernova mechanism. 

In the following, we briefly summarize the current status of our MRI project, in which 
we perform a local shearing box simulation to study the MRI-driven turbulence and their 
nonlinear properties bearing in mind the application to the supernova core (Masada, Taki- 
waki, and Kotake in preparation). The final goal of this project is to determine how much 
the conversion from the rotational energyf^l to the thermal ener g y occ urs via the MRI tur- 
bulence, which was treated parametrically in iThompson et al.l ( 120051 ). As will be shown 
later, our preliminary results suggest that the additional energy supply would affect the 
neutrino-driven explosion in the case of rapidly rotating cores at late postbounce phase. 

To study the nonlinear properties of the MRI, we solve viscous M HD equations by a 
finite-difference code that was originally developed by Sa no et al. ( 19981) . The hydrodynamic 



part is based on the second-order Godunov scheme (Ivan Leerl 119791 ). which consists of La- 
grangian and remap steps. The exact Riemann solver is modified to account for the effect 
of tange ntial magnet ic fields. The field evolution is calculated with Consistent MoC-CT 



method (jClarkd Il996l ) . The energy equation is solved in the conservative form and the vis- 
cous terms are consistently calculated in the Lagrangian step. The advantages of our scheme 
are its robustne ss for strong shocks and the satisfaction of the divergence-free constraint of 
magnetic fields (jEvans fc Hawleylll988l ; IStone fc Normanlll992l ). 



Rigorously speaking, 3D simulations are absolutely required for the study of the MRI 
because only by them the disruption of channel structures, mo de couplings, and the growt h 



of parasitic instabi li ties c an be accurately determined (e.g., iPessah &: Goodman! ( 120091 ); 
Longaretti Lesurl ( I20101 )). But in the following, we are only able to show our 2D results 
assuming axisymmetry, which are currently being updated to 3D. Bearing these caveats in 
mind, we briefly illustrate fundamental properties of the MR@ and discuss their possible 
impacts on the explosion dynamics. 



29 (that the differential rotation taps) 

30 which is not always familiar with supernova modellers, 
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t/t Iot k = k/(2n/L) 

Fig. 19. — (a) Temporal evolution of the volume averaged magnetic energy, Maxwell and 
Reynolds stresses in the fiducial run are represented by thick, dashed and dash-dotted curves 
respectively. The vertical axis is normalized by the initial magnetic energy -Z? mag ,o = Bq/8tt, 
that is (E mag ) = (E nmg ) / E maugfi , (M r(j> ) = (M r<t> ) / -E mag ,o, and (R r(j> ) = (R r(j) ) / E m3i g,o where 
single bracket denotes the volume average of the physical variables. Horizontal axis is normal- 
ized by the rotational period t rot = 2n/Q. Three typical evolutionary stages, (i) exponential 
growth stage, (ii) transition stage, and (iii) nonlinear turbulent stage, are denoted by white, 
dark gray and light gray shaded region, (b). The amplitude of power spectra of the magnetic 
energy \B(k)\ 2 along the k x and k z axes. The dashed and dotted curves describes the k x and 
k z components respectively. The vertical axis is normalized by the initial magnetic energy 
Em&g,o- The wave-numbers shown in the horizontal axis are normalized by the 2n/L. The 
upper thick curve demonstrates the Kolmogorov slope for isotropic incompressible turbu- 
lence oc At 11 / 3 . The lower slope is proportional to k~ 6 just for the reference. In our local 
simulations, we prepare a 2D shearing box which is threaded by the initial vertical field 
B = 1.4 x 10 12 G with imposing initial gas pressure P = 4 x 10 28 dyn cm -2 , angular veloc- 
ity Qq = 2000 sec -1 and shear parameter q = 1.25, respectively. Here the shear parameter 
represents the degree of differential rotation as q = — dlnQ/dlnr with r being the radial 
coordinate. The other parameters are fixed to mimic the properties in the vicinity of the 
PNS's surface in the postbounce phase (i.e. the density is taken to be po — 10 12 g cm -3 and 
the box size is taken to be L x = 4 km (horizontal) x L z = 1 km (vertical)). 
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Fig. 20. — Four snapshots representing the temporal evolution of MRI from the linear to 
the non-linear phase. Panels (a), (b), (c) and (d) correspond to the time slice t = 3t rot , 
4£ rot , 6t TOt and 20t rot respectively. The color bar indicates the size of the toroidal magnetic 
field normalized by the strength of the initial vertical magnetic field (B = 1.4 x 10 12 G). 
Note that the vertical and horizontal are normalized by the typical spatial scale of the local 
portion in the proto-neutron stars which is selected as L = 1km in the fiducial run. Three 
typical evolutionary stages are observed, (i) exponential growth stage, (ii) transition stage, 
and (iii) nonlinear turbulent stage (compare with Figure [TH]). 
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Fig. 21. — The dependence of the time- and volume- averaged (a) Maxwell (green squares) 
and Reynolds stresses (orange diamonds) and (b) magnetic energy (blue circles) and in the 
turbulent state on the strength of the initial magnetic energy E ini . Both the vertical and 
horizontal axes are normalized by the initial magnetic energy of the fiducial model E maig ^. 
The initial field strength is varied from ~ 10 11 to ~ 6 x 10 12 G. The dashed lines denote 

1 /2 

reference slopes proportional to (a) E„[ ag and (b) E ma _g respectively. We take the time-average 
of volume-averaged quantities during 90t rot < t < 100t rot . 



3.0.4- Local simulations in axisymmetry 



As it is well known, there are three typical evolutionary stages observed in our local 
simulations. Th ey are (i) exponential growth stage, (ii) tra nsition stage, and (iii) nonlinear 
turbulent stage (IHawley fc Balbudll992l ; lHawley et al.lll995l ). Each stage is denoted by white, 
dark gray and light gray shaded regions in Figure | 



In the stage (i), the channel structure of the magnetic field exponentially evolves and 
inversely cascades to the larger spatial scale with small structures merging as the magnetic 
field is amplified. This is beca use the channel modes of the MRI are a exact solution for 
the nonlinear MHD equations (IGoodman &: Xul Il994l ) . The temporal evolution of channel 
structures for the toroidal magnetic field and their inversely cascading nature are shown 
in panels (a) and (b) of Figure [20j Then in the stage (ii), the channel structure of the 
magnetic field is disrupted via the parasitic instability and mag netic reconnection at the 
transition stage (IGoodman fc Xulll994J : lObergaulinger et al.ll2009l ). The channel disruption 
induces a drastic phase-shift from the coherent structure to the turbulent tangled structure 
of the field as is illustrated by panel (c) of Figure [20j The magnetic energy stored in the 
amplified magnetic field is then converted to the thermal energy of the system. Finally in 
the stage (iii), the nonlinear stage emerges after the channel structures are disrupted to be a 
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turbulence state driven by the MRI (see panel (d) in Figure [20]). Figure [T9k shows that, at 
this stage, the turbulent Maxwell stress [M r ^ = —B r B^/A7r) dominates over the Reynolds 
stress {R rt j> = pVrdv^). Therefore the turbulent Maxwell stress is the main player to convert 
the free energy stored in the differential rotation into the thermal energy of the system. 

Figure [T9b shows the power spectra of the magnetic energy in the nonlinear turbulent 
stage. The logarithmic slope of the spectrum roughly coincides with —11/3 in the regime 
10 ^ ^ ~ 100 while it steepens at higher wave-numbers where numerical effects become 
important (like k~ & ). This indicates that the turbulent state in the MRI can be represented 
by the Kolmogorov spectrum for a homogeneous, incompressible turbulence that scales as 
fc-n/3 These features wo uld be consistent with the 3D properties of the MRI-driven turbu- 
lence (IHawley et al.lll995l ). It has been pointed out in the previous axisymmetric local shear- 
ing box simulation that the channel structures d o neither disrupt nor transit to the turbulent 



state without a relatively strong resistive effect (ISano fc Inutsukall200ll ; lObergaulinger et al. 



20091 ). We however adopted the simulation box with a large aspect ratio (L x /L z = 4) and 
quasi-isothermal e quation of state O y ~ 1), both of which could facilitate to disrupt the 
channel structures (IBodo et al.ll2008l ) as in 3D simulations. 



Figure [21] shows the time- and volume-averaged Maxwell (green squares) and Reynolds 
stresses (orange diamonds, panel (a)), and magnetic energy (blue circles, panel (b)) in the 
turbulent stage as a function of the initial magnetic energy. It is shown that the Maxwell 
stress depends on the initial magnetic flux and it provides the following scaling relation, 



((M r4 ,)) ex E 



mag, inn 



(4) 



with -Emag^ni = Bf ni /8n. This suggests that the magnetic flux initially penetrating the 
system controls the nonlinear transport properties of the MRI-driven turbulence. This is 



qualitatively con s istent with the previous res ults in the accretion disk (IHawley et al.lll995 



Sano et al.l Il998l . l2004j ; iPessah fc Chanl 120081 ) while the power-law index obtained here is 
slightly larger. Note that the magnetic energy at the nonlinear turbulent stage has a linear 

1 /2 _____ 

dependence on the initial magnetic flux (because ({E mae )) oc P magi i m t ^ Ani (Figure l2~TT b))). 
This is also qualitatively consistent with Hawley et al. (1995). In this way, we investigated 
parameter dependence of the rotational shear (q parameter) as well as the initial pressure 
(P) on the saturated value of the Maxwell stress and deduced a scaling relation that yields 
M r< p oc B 2 P l / A [q/ (2— g)] 1//4 f2 1//2 . Here the shear parameter represents the degree of differential 
rotation as q = — dhxQ/dlnr with r being the radial coordinate. 

Using the scaling relation obtained in our local simulations, the heating rate maintained 
by the MRI turbulence can be estimated as, 



e M Ri = M l4> qQ 
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= /( J B pb / J B ) 2 (Ppb/ J Po) 1/4 4 / b 4 ( fi Pb/^o) 1/2 ] 
~ 10 30 erg cm -3 sec -1 x 

, 1/4 ( B ph \ 2 f Ppb / ^pb\ 



(5) 



3/2 



where /30 represents a typical amplification factor of the stress normalized by 30 (e.g., Figure 
EE]), <Zpb,i is the shear rate normalized by unity and s p b,i is the ratio of the vorticity to 
the shear normalized by unity. Since we are interested in the growth of the MRI in the 
postbounce phase, the original magnetic flux, gas pressure angular velocity, and the shear 
rate B = 1.4 x 10 12 G, P = 4 x 10 28 erg cm" 3 , Q = 2000 rad s _1 , and q = 1.25 are replaced 
by Bph, Ppb, ^ p b and g p b with some amplification factors. 

A typical volume in which the MRI-driven turbulence is active, may be estimated as 
Vmbi = 4:irR 2 h ~ 10 21 cm 3 P 2 /i 6 , where R 7 = R/10 7 cm is typical radius of the PNS nor- 
malized by 10 7 cm, and = /z/10 6 cm is the typical radial width of the MRI-active layer 
normalized by 10 6 cm. Then the energy releasing rate Lmri = cmriHiri becomes, 



>MRI 



10 51 P^ 6 erg sec" 1 



x 



(6) 



t 1/4 

J30 9pb,lS pb ,l 



B 



P b 



(■ 



P 



P b 



iop y ViooPq 



1/4 



n 



P b 



n j 



s 3/2 



Typical timescales of neutrino-d riven explosions obser v ed in recent super nova simulations 
are > 400 ms after bounce (e.g., iMarek fc Jankal (120091 ); ISuwa et al.l fj2010h ) . So the energy 
deposition could be as high as 10 50 erg if the core rotates rather rapidly as taken in the above 
estimation. If this is the case, the turbulent heating due to the MRI is expected to play an 
important role for assisting the neutrino-driven explosion. 

It should be noted that nonaxisymmetric mo des of the parasitic instabil i ty could disrupt 
the ch annel structures more effectively than in 2D ( Pessah fc Goodmanll2009 ; Longaretti fc Lesur 
20101 ). To obtain more accurate estimate of the amplification factor of / (equation (jSJ)), 3D 
simulations are unquestionably required, which we are currently undertaking. 



3.1. Gravitational waves 



As already mentioned in section 12.21 rapid rotation, necessary for producing strong 
bounc e GW signals, is likely to obtain ~ 1% of massive star population (e.g- JWoosley Bloom 
( 120061 )). This can be really the case, albeit minor, for progenitors of rapidly rotating metal- 
poor stars, which experience the so-called chemically homogeneous evolution ( IWoosley fc Heger 
20061 ; lYoon &; Langerl l2005bl ) . In such a case, the MHD mechanism could work to produce 
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energetic explosions, which are receiving great attention rece n tly as a possible r e levan ce to 
magnetars an d collap s ars fe.g.. iMacFadven &: Wooslevl ( 119991 ); iMacFadyen et al.l ( 120011 ) . see 
Dessart et al.l ( 120081 ); lHarikae et al.l (120091 |2010| ) for collective references), whic h are pre- 



sumab ly linked to the formation of long-duration gamma-ray bursts (GRBs) 
( 120061 ) for review). 



e.g. 



Meszaros 



Among the previous studies focusing on the bounce signals (e.g., references in section 



V 1 " 

MHD mechanism ( 


Yamada & Sawai 


20041; 


Kotake et al. 


2004b 


Obergaulinger et al. 


2006b: 


Cerda-Duran et al. 


2007 


Shibata et al. 


2006|; 


Scheidegger 


et al. 


2010a 


). This may be because 



the MHD effects on the dynamics as well as their influence over the GW signals can be 
visible only for cores with p recollapse magnetic fields over B > 10 12 G (lObergaulinger et al. 



2006b 



Kotake et al.l l2004bl ) . Conside ring that the typical ma gnetic-field strength of GRB 
progenitors is at most ~ 10 11-12 G (jWoosley fc Heger 2006 ). this is already an extreme 
situation. In a more ext remely case of Bo ~ 10 13 G, a secularly growing feature in the 



wavefo rms was observed (lObergaulinger et al.ll2006al ; IShibata et al.ll2006l ; IScheidegger et al. 
2010ai ). In the following, we summ arize the GW signatures based on our 2D SRMHD 



simulations ( jTakiwaki fc Kotakell201ll ) and discuss how a peculiarity of the MHD mechanism 
could be imprinted in the GW signals. 



The left panel of Figure [22] shows the gravitational waveform with the increasing trend, 
which is obtained for a model with strong precollapse magnetic field (B = 10 12 G) also with 
rapid rotation initially imposed (f3 parameter = 0.1 %). Such a feature cannot be observed 
for a weakly magnetized model (Bq = 10 11 G. right panel). To understand the origin of the 
increasing trend, it is straightforward to look into the quadrupole GW formula, which can 
be expressed as 



hJ7(x,t) e=2 ^ =0 Ia* 3 

R 



1 20 



i rrE2,20, 

where 1^ ' { 



is 



rjiLL 



£2,20/ 



15 • 2n 

— sin 0, 




7T 



[e.g., 



Thornel (Il980h ). Af 2 can be expressed as 



A E2 
^20 



A 



E2 



+ A 



E2 



+ A 



E2 



(7) 



(8) 



(9) 



1 20 (hyd) ~ ^20 (mag) ~ ^20 (grav) 

On the right hand side, the first term is related to anisotropic kinetic energies which we call 
as hydro dynamic part, 

r dr J20 (hyd) 



4 E2 

^20 (hyd) 



rE2 

J 20 (hyd) 



da 

< x/15 Jo 

r 2/„. 2/o,,2 



p,W 2 (v r 2 (3^ 2 - 1) + v e 2 {2 - 3/x 2 ) - - Qv r v e 



(10) 
(11) 
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the second term is related to anisotropy in gravitational potentials which we call as the 
gravitational part, 

G32vr 3 / 2 f 1 . r°° 



aE2 

^20 (grav) 



rE2 

J 20 (grav) 



(12) 



X 



(13) 



-rd r $(3// 2 - 1) + 3d e $ fiy/l - p 2 , 

and finally the third term is related to anisotropy in magnetic energies that we refer to as 
magnetic part, 

dp I r 2 rfr/ 2 E 2 (mag) , (14) 

(15) 



A E2 

^20 (mag) 



J 



fE2 



20 (mag) 



[6 r 2 (3/i 2 - 1) + 6 e 2 (2 - 3^ 2 ) - 6/ - Qbrhp^l - p 2 } 
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Fig. 22. — Gravitational waveforms with an increasing trend (left) or not (right panel). 
Initial rotation parameter is both set to be /3 = 0.1 %, while the precolla pse magnetic field is 
taken as 10 12 G (left panel) and 10 11 G (right panel), respectively (from iTakiwaki fc Kotake 
( 1201 ll )). The total wave amplitudes are shown by the red line, while the contribution from 
the magnetic fields and from the sum of hydrodynamic and gravitational parts are shown by 
blue and green lines, respectively (e.g., equation (ITS"]) and equations (lll|13p ). Note that the 
bounce GW signals ((£ — £& < 20 ms) are not affected by the magnetic fields significantly, 
and they are categorized into the so-called type I or II waveforms. 



The right panel of Figure [23] shows contributions to the total GW amplitudes (equa- 
tion (6)) for the strongly magnetized model (left panel of Figure [221 . in which the left- 
hand-side panels are for the sum of the hydrodynamic and gravitational part, namely 
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Fig. 23. — Left panel shows distributions of entropy [/c^/baryon] (left) and logarithm of 
plasma (3 (right) for model B12X1/30.1 at 100ms after bounce. The white arrows in the left- 
hand side show the velocity fields, which are normalized by the scale in the middle left edge 
(0.5c). Right panel shows the sum of the hydrodynamic and gravitational parts (indicated 
by "hyd and grav" in the left-hand side) and the magnetic part (indicated by "mag" in the 
right-hand side), respectively. The top and bottom panels represent the positive and negative 
contribution (indicated by (+) or (-)) to Afg, respectively (see text fo r more detail). The 
side le ngth of each plot is 4000(km)x8000(km). This figure is taken from lTakiwaki fc Kotake 

feonh . 
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log ± 



xE2 
J ' 



20 (hyd) 



4- f E2 

' J 20 (grav) 



(left top(+)/bottom(— ) (equations (jll|13p ). and the right-hand- 
side panels are for the magnetic part, namely log ^±/20 2 (mag)) (right top(+)/bottom(— )) 
(e.g., equation (fl5l)). By comparing the top two panels in Figure 19, it can be seen that 
the positive contribution is overlapped with the regions where the MHD outflows exist. The 
major positive contribution is from the kinetic term of the MHD outflows with large radial 
velocities (e.g., +p 1r W 2 v r 2 in equation ( FIT]) ). The magnetic part also contributes to the pos- 
itive trend (see top right-half in the right panel (labeled by mag(+))). This comes from the 
toroidal magnetic fields (e.g., +b ( p 2 in equation (lisp), which dominantly contribute to drive 
MHD explosions. 

Figure [24] shows the GW spectra for a pair models of B12X5/31 and B11X5/31 like 
in Figure [22] that with or without the increasing trend. Regardless of the presence (left 
panel) or absence (right panel) of the increasing trend, the peak amplitudes in the spectra 
are around 1 kHz. This is because they come from the GWs near at bounce when the 
MHD effects are minor. On the other hand, the spectra for lower frequency domains (below 
~ 100 Hz) are much larger for the model with the increasing trend (left panel) than without 
(right panel). This reflects a slower temporal variation of the secular drift inherent to the 
increase-type waveforms (e.g., Figure I2"2]) . It is true that the GWs in the low frequency 
domains mentioned above are relatively difficult to detect due to seismic noises, but a recently 



(Kawamura 


2006; 


Kudoh et al. 


2006) 



in Figure |24"|) . Our results suggest that these low-frequency signals, if observed, could be one 
important messenger of the increase-type waveforms that are likely to be associated with 
MHD explosions exceeding 10 51 erg. 
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Fig. 24.— Gravitational-wave spectrum for representative mod els with the ex- 
pected detection limit s of the first LIGO (lAbbott et al.l I2005I ). the advanced 



LIGO (|Weinsteinl I20021). Large-scale Cryogenic Gravitational wave T elescope (LCGT) 



Kudoh et al 



(IKuroda fc LCGT Collaboration! l2010f ). and Fabry-Perot type DECIGO (IKawamural I2006 



2006) . It is noted tha t /i c h ar is the characteristic gravitational wave strain de- 



fined in lFlanagan fc Hughes! (11998! ). The supernova is as s umed to be located at the distance 
of 10 kpc. This figure is taken from iTakiwaki fc Kotakd (I201ll ). 
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3.2. Neutrino Signals 



As of now, SN1987A remains the only astrophysical neutrino source outside of our solar 
system. Even though the detected events were just two dozen, these events have been studied 
extensively (yielding ~ 500 papers) and have allo wed us to have a confid ence that ou r basic 



pictu re of the supernova physics is correct (e.g., iSato Suzuki! (Il987l ). see iRaffeltl (12012 
20021 ) for a recent review). For a next nearby event, it is almost certain that the flagship 
detectors like Super-Kamiokande and IceCube are able to detect a SN "neutrino light curve" 
with high statistics. For the detectors, the horizon to the sources now extends out to about 
100 kpc and thus covers our galaxy and its satellites. A future megaton-class detector reaches 
as far as Andromeda galaxy at a distance of 780 kpc (e.g., Hyper-Kam iokande, Memphys, 
and LBNE) and large-scale scintillator (e.g., HALO (E ngel et al.l 120031) that is an upgrade 



of SNO) or liquid- Argon detectors (GLACIER ( Autiero et al.ll2007Q that is an upgrade of 
ICARUS) will also play an important role (see Scholberg ( 201ol )) for a complete list and 
details). 

Before core bounce, neutrinos of different flavors are initially tapped in their relative 
neutrino spheres with increasing their Fermi energies as the central density becomes higher. 
After bounce, the shock dissociates nuclei into free nucleon, which drastically increases the 
number of protons, leading to a sharp increase in electron capture and production of the 
prompt u e burst. At this neutronization epoch (in the first 10-20 ms postbounce), the largest 
difference i n the light-curves amon g v e and the rest of the neutrino species can be seen (e.g., 
Figure 1 in lKachelriefi et al.l (120051 ) ). Since neutrino oscillations take pla ce only when there is 
a difference in the neutrino fluxes among different species (see section 3 iKotake et al.l (120061 ) 
for more details on the neutrino oscillations), the oscillation effects in the neutronization 
epoch would be strongest. Unfortunately, however, they are very hard to measure because 
the currently-running detectors are primarily sensitive to v e signals that are produced by the 
inverse beta process v e + p — > e + + n. It should be noted that the detection the v e bursts is 
important for studying not only neutrino oscillatio ns but also for determining the direction 
and distance to the source ( KachelrieB et al. 20051). In a megaton-class C herenkov detector 
in which gadolinium is added to catch neutrons (IBeacom &: Vagina 120041 ). and large liquid- 
Argon detectors (like ICARUS and GLACIER cited above) are expected to be powerful v e 
detectors. 

During the subsequent accretion phase (approximately before 1 s postbounce), the 
bounce shock stagnates and matter falls in through the stalled shock to the center, releasing 



gravitational energy that powers neutrino e mission. For r e 



in the mass range of 1 0.8 to 25 ( e - P-- iKitaura et al.l (12 006) ; 



KachelrieB et al.l (l2005l ): lHudepohl et al.l (l2010h : 



ative ly well-studied pro g enitor s 



Fischer et al 



Thompson et al.l (120031 ); 



(120 lOl ) and references therein) , 
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one generally finds L Ve ~ L Pe > L Vx and (E Ux ) > (E Pe ) > (E Ue ) with L Vi and (E Vi ) represent- 
ing luminosity and average neutrino energy for the neutrino species of i — v P , u P , v x ^v^v T 



and their anti-particles). The accretion phase is followed by the PNScj cooling in which the 
accretion stops and the PNS settles to become a neutron star after a successful revival of the 
stalled bounce shock into explosions. Since the explosion mechanism is still under debate 
(see discussion in section [2$, it is at present difficult to tell the accurate turn-over time. 
But it sho uld be less than 1-2 s po stbounce, otherwise the central PNS would collapse into a 
black hole (jO'Connor fc Ottll201ll ). In the PNS cooling phase, the luminosities of all neutrino 
species become similar L (with the energy hierarchy of (E Ux ) > (E De ) > (E Ve )) 

and decrease monotonically with time. Due to the similarity in the luminosities and spectra 
among v e and u x and also due to the darkening with time, it is much harder to see flavor 
oscillation effects by the currently running detectors. 

Neutrinos streaming out of t he iron core interact with m atter via the Mikheyev-Smirnov- 
Wolfenstein (MSW) effect (e.g.. iMikheev fc Smirnov! (119861 )) firstly in propagating through 
progenitor envelope and then through the Earth before reaching detectors. Such effects 
have been extensively investigated so far from various points of v iew, with a focu s such a s 



on the progenito r 



dep endence of the early neutrino burst (e .g., 



KachelrieB et al.l (20051 )) and on the Earth matter effects (e.g.. iLunardini &: Smirnov! (120011 ); 



Takahashi et al 



(20QM); 



Dighe et al.l (120041 ) ). The conversion efficiency via the MSW effect depends sharply on den- 



sity, equivalently electron fraction gradients, thus sensitive to discontinuities produced by 



Schirato & Fuller! ( 


2002); 


Lunardini & Smirnov 


(2003h; 


Tomas et al. 


(2004 


); 


Daseupta & Diehe 


(2007 


); 


Lunardini & Peres! ( 


2008l).see 


Duan & Kneller 


'2009h for a recent review"). The shock 



passage to the so-called high-resonance regions may be observed as a sudden decrease in v e 
events in the case of normal mass hierarchy (or v e in the case of inverted mass hierarchy). 
Such features monitoring time evolution of the density profile like a tomo graphy, thus could 



provid e a powerful test o f the m i xing angle and the mass hierarchy (e.g . , iDighe fc Smirnov 



(120001 ): iTakahashi et all (I2003bf ): iTomas et all (120041 ): iFogli et al.l (120051 ) ). 



It is rather recently that the importance of collective flavor oscillations was widely 
recognized. Neutrinos streaming out of the neutrino spheres are so dense that they provide 
a large matter effect for each other. The collective effect usually takes place between the 
neutrino sphere and the mentioned MSW region. Regardless of the inherent non-linearity 
and the presence of multi-angle effects, the final outcome for the emergent neutrino flux 
seems to be converging after a series of extensive study over the past years at least for the 



31 (proto-neutron star) 
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ID models. Alth ough the collective e ffects will not be significant to assist the neu trino-driven 



explos i ons fe.g.. IChakrabortv et al.l ( 1201 ll ) ; iDasgupta et al.l ( 1201 ll ) , see however ISuwa et al. 



( 120 111 ); iPejcha et al.l ( 120111 ) ), they are predicted to emerge as a distinct observable feature 



in their energy spectra (see iRaffelt fc Smirnovl 120071 ; iDuan et al.ll2008l . |201Q| ; lDasguptall2010 



for reviews of the rapidly growing research field and collective references therein). 

An important lesson from SN1987A is that for explaining the duration of the events, 
there was not other energy-loss channel b ut for the ordinary n eutrinos in the context of the 
Standard Model of particle physics (e.g., iRaffeltl (119961 120021 )). The next SN event could 
provide an opportunity to study also a nontrivial property of neutrinos, such as the magnetic 
dipol e moments. The re s onant s pin-flavor convers i on ha s been also studied both a nalyt ically 



Schechter fc Valid (119811 ): iLim fc Marcianol (119 88); 



Akhmedov fc Khlopovl (ll988f ) and 



(e.g., 

references therein) and numerically (e.g., lAndo &; Satol ( 120031 )). which can transform some 



of the prompt v e burst into v e in highly magnetized supernova envelopes, leading to a huge 
v e burst. The mentioned important ingredients related to the flavor conversions due to the 
MSW effects, collective effects, and the electromagnetic effects in the supernova environment 
have been studied often one by one in each study without putting all the effects together, 
possibly in order to highlight the new ingredient. In this sense, all the studies mentioned 
above should be regarded as complimentary towards the precise predictions of SN neutrinos. 

Here it should be mentioned that most of those rich phenome nology of SN neutr i- 



nos have been based on ID spherically symmetric si mulations (e.g., Totani et al.l ( 1998) 



( 2005 ): 

(bonh : 



Lunardini et a 



Lunardini Sz Smirnov J 20031); Takahashi et al. J2003ah : |^mas et al.l (EoO^ 



JgOOgJ); IDuan et al.l f|2008h : Irliidepohl et~all fj201oh : Icherrv et al" 



Kachelriefi et al. 



Sarikas et al.l (120111 ) and references therein). Apar t from a simple parametriza - 



tion to mimic anisotro p y and Stochasticit y of the shock (e.g.. iFriedland fc Gruzinovl (120061 ): 
Choubey et al.l (120071 ): iReid et al.l (120111 )). there have been not so many studies focusing 
how global asymmetry of the explosion dynamics obtained in recent super nova simulations 
(e.g., s ection [2] and El) would affect the neutrino oscillations (see, however, IDuan fc Kneller 
(120091 ): IDasgupta et al.l (120081 ) ). This is mainly due to the lack of multi-D supernova models, 
which are very computationally expensive to continue the simulations till the shock waves 
propagate outward until they affect neutrin o transformations. It is only recently t hat several 
studies along this line have been reported (ILund et al.ll2010l : iBrandt et al.l 1201 if ), in which 
the collective or MSW effects are rendered t o be treated in an app roximate manner. By an- 
alyzing the 2D results of a 15 M model by iMarek fc Jankal (120091) who incl u ded o ne of the 
best available neutrino transfer approximations (e.g., Table 1). ILund et al. ( 2010 ) pointed 
out that fast time variations caused by convection and SASI lead to significant modulations 
around a few hundred Hz, which can be visible in IceCube or future megat on-class de t ectors 
for the galactic SN source. Based on the 2D results of 20 M & models by lOtt et al.l (120081 ) 
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who reported the first 2D multi-angle transport simulations, iBrandt et al.l ( 1201 ll ) also ob- 
tained the similar results. From their rapidly rotating model, they also pointed out that 
a rapid rotation of p recollapse SN cores imprints strong asymmetries in the neutrino flux 



dKotake et al.ll2003al) as well as in its light curves (IJanka &; Monchmeyerlll989t IBrandt et al. 
201lh . 



In the following, we briefly summarize our findings (IKawagoe et al.ll2009l ) in which we 
studied exploratory the neutrino oscillations in the context of the MHD mechanism. As 
mentioned in section [31 it is recently possible for special relativistic simulations to follow the 
dynamics of MHD explosions continuously, starting from the onset of gravitational collapse, 
through core-bounce, the ma gnetic shock-reyival, till the shock-propagation to the stellar 
surface. Based on our models ( ITakiwaki et al.ll2009l ). we calculated numerically the neutrino 
flavor conversion in the highly non-spherical envelope throug h the MSW effect. Th e neu- 
trino transport was simply treated by a leakage scheme (e.g., ITakiwaki et al.l (120091 ) ) . The 
emergent neutrino spectra was assumed to take a Fermi-Dirac type distribution function 
and the neutrino temperature was estimated to take the average matter temperature on the 
neutrino sphere. As explained section 3.0.2, the density profile along the polar and equato- 
rial direction is very different due to the strong explosion anisotropy inherent to the MHD 
explosions. In the following, we pay attention to the anisotropic shock effect on the MSW 
effect. As will be discussed, we could observe a sharp dip in the neutrino event only seen 
from a polar direction, albeit depending on the mass hierarchy and the the mixing angle 
of 6* 13 . An advantage of the MHD models is that the shock revival can occur much faster 
after bounce compared to the other proposed mechanism. Therefore the neutrino luminosity 
could remain higher than those in the other mechanisms, which could potentially enhance 
the detectability due to the early shock arrival to the resonance region. For simplicity, the 
effects of neutrino self-interactions are treated very phenomenologically and the resonant 
spin-flavor is not considered. Even though far from comprehensive in this respect, we pre- 
sented the first discussion how the magneto-driven explosion anisotropy has impacts on the 
emergent neutrino spectr a and the result i ng ey ent number observed by the SK for a future 



Galactic supernova (e.g., Kawagoe et al.l ( 120091 )). 



3.2.1. Possible neutrino signatures from MHD explosions 

One prominent feature of the MHD models is a high degree of the explosion asphericity. 
Figure [25] shows several snapshots featuring typical hydrodynamics of the model, from near 
core-bounce (0.8 s, left), during the shock-propagation (2.0 s, middle), till near the shock 
break-out from the star (4.3 s, right), in which time is measured from the epoch of core 
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Fig. 25. — Snapshots showing MH D explosions of CCSN e at 0.8, 2.0 and 4.3 sec after 



core bounce from left to right (e.g., iTakiwaki et al.l (120091 ) ) . In each figure, contour of the 
logarithmic density [g/cm 3 ] (left) and entropy per baryon [&b] (right) are shown. The unit 
of the horizontal and the vertical axis is in [cm] . Note that the difference of the length scale 
for each panel and that the outermost radius of the progenitor is ~ 3 x 10 10 [cm]. 

bounce. It can be seen that the strong shock propagates outwards with time along the 
rotational axis. On the other hand, the density profile hardly changes in the equatorial 
direction, which is a generic feature of the MHD explosions. 

As well known, the flavor conversion through the pure-matter MSW effect occurs in the 
resonance layer, where the density is 

, , 3 / Am 2 \ /10MeV\ /0.5\ „ , s 

p re , ~ 1.4 x 10 W ( TF5 ^j (— ) (^j cos20 (16) 

where Am 2 is the mass squared difference, E u is the neutrino energy, Y e is the number of 
electrons per baryon, and 9 is the mixing angle. Since the inner supernova core is too dense 
to allow MSW resonance conversion, we focus on two resonance points in the outer supernova 
envelope. One that occurs at higher density is called the H-resonance, and the other, which 
occurs at lower density, is called the L-resonance. Am 2 and 9 correspond to Am 2 3 and #13 
at the H-resonance and to Am 2 2 and 9\i at the L-resonance. 

Figure [26] are evolutions of the density profiles in the polar direction every 0.4 s as a 
function of radius. Above and below horizontal lines show approximately density of the 
resonance for different neutrino energies which are 5 and 60 MeV, respectively. Blue lines 
(left panel) show the range of the density of the H-resonance, and sky-blue lines (right 
panel) show that of the L-resonance. Along the polar axis, the shock wave reaches to the 
H-resonance, ~ O(10 3 )g/cm 3 at ~ 0.5 s, and the L-resonance, ~ 0(l)g/cm 3 at ~ 1.2 s. It 
should be noted that those timescales are very early in comparison with the ones predicted 
in the neutrino-driven explos ion models , typic a lly ~ 5 s and ~ 1 5 s for the H- and L- 



resonances, respectively (e.g., iFogli et al.l ( 120051 ): iTomas et al.l ( 120041 )). This arises from the 
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fact that the MHD explosion is triggered promptly after core bounc e , whi c h is in sharp 



contra st to the neutrino-driven delayed explosion models ( IFogli et al.l (120051 ) ; iTomas et al. 



( 120041 ) ). The progenitor of the MHD models, possibly linked to lo ng-duration gamma-ray 
bursts, is more compac t due to a chemically homogeneous evolution (jYoon fc Langerll2005bl ; 
Woosley fc Hegerll2006l ). which is also the reason for the early shock-arrival to the resonance 
regions. 
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Fig. 26. — The density profile in the polar direction every 0.4 s as a function of radius, (a), 
(b) and (c) correspond to 0.4, 0.8 and 2.8 s, respectively. The horizontal blue lines show 
approximately density of the H-resonance for different neutrino energies which are 5 (above) 
and 60 MeV (below), respectively. The green line is an original density profile for 2.0 s. Same 
as the left panel, but for 1.2 (d), 2.0(e), and 3.2 (f) s, respectively. Note that a sharpness of 
the shock along the polar direction is modified from the one in the MHD simulations (green 
line showing an original density profile for 2.0 s), which is made to be inevitably blunt due 
to the employed nume rical scheme using an artificial viscosity to capture shocks (taken from 
Kawagoe et aD ( hood )). 



Figure [27] is the energy spectra of v e for (a), (b) and (c) in the case of inverted mass 
hierarchy. The solid lines (red lines) and the dotted lines (blue lines) are the spectra of the 
polar direction and the equatorial direction, respectively. At the panel (a), an enhancement 
of the low-energy side of the polar direction is seen, which is as a result of the supernova 
shock reaching to the H-resonance region. The survival probability of v e can remain non- 
zero when the ste ep decline of th e dens ity at the shock front changes the resonance into 
non-adiabatic (see iKawagoe et al.l (120091 ) for more details). As the shock propagates from 
the high density region to low energy region, the shock effect transits from low-energy side to 
high-energy side in spectra, because the density at the resonance point, p res , is proportional 
to E~ x . In fact, the energy spectrum of the polar direction becomes softer in the high-energy 
side than for the equatorial direction as shown in Figure |27J (b). 
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Figure [28] shows the time evolution of the event number in SK, where the solid line 
(red line) and the dotted line (blue line) are for the polar and the equatorial direction, 
respectively. The shock effect is clearly seen. The event number of the polar direction shows 
a steep decrease, marking the shock passage to the H-resonance layer. The change of the 
event number by the shock passage is about 36% of the event number without the shock. 
Since the expected events are ~ 2500 at the sudden decrease (Figure 128]) . it seems to be 
quite possible to identify such a feature by the SK class detectors. Such a large number 
imprinting the shock effect, possibly up to two-orders-of magnitudes l arger than the ones 
predicted in the neutrino-driven explosions (e.g., Fig 11 in lTomas et al.l (12004J )). is thanks to 
the mentioned early shock-arrival to the resonance layer, peculiar for the MHD explosions. 
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Fig. 27. — v e spectra at the surface of the star in the case of (a), (b) and (c). We assume the 
inverted mass hierarchy and sin 2 2# 13 = 10 -3 . The solid lines (red lines) and the dotted lines 
(blue lines) are the spectra of polar direction and equatorial direction, respectively. Note in 
this figure that the neutrino luminosity and spectra are modeled to be the same as those 
obtai ned in a ID full-scale numeric al simulation by the L awrence Livermore group (jWilson 
1985I ). These figures are taken from lKawagoe et al.l (]2009l ). 
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Fig. 28. — Same as Figure 1271 but for the expected event number of u e (left panel) in the 
SK as a function of the time. The solid line (red line) and the the dotted line (blue line) 
is for the polar and the equatorial direction, respectively. The supernova is assumed to be 
located at the distance of 10 kpc. The error bars represent the la statistical errors only. 
Right panel is for the ratio of high- to low-energy events which can be a useful quantity to 
characterize th e chan ge (dip in the right panel) of neutrino signals by the shock effect (see 
Kawagoe et al.l (120091 1 for details). 



Finally we exploratory discuss possible impacts of the spectral swap between and u^ , 
which is one of a probable outcome of the collective neutrino oscillations (jFogli et al.l 120071 ). 
In our model, the original neutrino flux of v e is larger than that of v x for the low-energy side, 
which is vice vers a for the high-e n ergy side. Due to the spectral swap, the above inequality 



sign reverses (see iKawagoe et al.l ((20091) for more detail). As a result, the event number in 
the polar direction is expected to increase by the influence of the shock. To draw a robust 
conclusion to the self-interaction effects , one needs t o perf orm a more sophisticated analysis 



such as the multi-angle approach (e.g., iDuan et al.l (120101 )). which is a major undertaking. 
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4. Summary and Discussion 
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hanism 
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Excess for equator 
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s: black-hole forming 


Burst signals 
(bounce & BH formation) 


Disappearing signals 


No photon (?) 








MHD mechanism 


Burst & tail 

(rapid rotation + 
magnetic fields) 


-Anisotropy in SK 
events (MSW effect) 

■~ e bursts (RSF) 


■ r-process cites ? 
■Path to 
hypernovae ? 



Table 2: A summary illustrating the candidate supernova mechanisms (vertical direction in 
the table) and their expected features of the emergent multimessengers (horizontal direction). 
Regarding the photon messengers, possible signatures associated with nucleosynthesis are 
only partly presented (see section 1 for a number of important findings regarding the light- 
curve and spectra modeling). 

The aim of writing this article was to provide an overview of what we currently know 
about the explosion mechanism, neutrinos, gravitational waves, and explosive nucleosynthe- 
sis in CCSNe to bridge theory and observation through multimessenger astronomy. Not to 
inflate the volume, we had to limit ourselves to focus primarily on the findings based on our 
numerical studies, which are thus limited to GWs (section 12. 2p and explosive nucleosynthesis 
(section |2. 3j) in the neutrino-heating mechanism (section [2]), and to GWs (section 13.1 1) and 
neutrino signals (section 13 .2 j) in the MHD mechanism (section [3]). Apparently our current 
presentation is weak on the photon side mainly because there is a big gap to connect the out- 
comes obtained in the multi-D neutrino-radiation-hydrodynamic simulations covering only 
the very centre of the star, to the light-curve and spectra modeling that requires multi-D 
photon-radiation-hydrodynamic modeling (far) after the shock-breakout of the star. The 
gap is not only the matter of physical scale, but also the matter of numerical difficulty 
Note that table [2] is not intended to show an overview covering everything concerning the 
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SN multi-messengers, which is far beyond the scope of this review. We boldly present the 
incomplete table here, hoping that the shortage might give momentum to theorists for fixing 
the problems and for gaining much more comprehensive multimessenger perspectives for the 
future. 

In the neutrino-heating mechanism, stellar rotation holds a key importance to char- 
acterize the SN multimessengers. The whole story may be summarized as follows. If a 
precollapse ir on core has a "canonical " rotation rate as p redicted by recent stellar evolution 
calculations (IMaeder fc Meynetl l2000l ; Irleger et al.l 120051 ) . collapse-dynamics before bounce 



proceeds spherically and the postbounce structures interior to the PNS are essentially spher- 
ical. Typically later than ~ 100 ms after bounce, the bounce shock turns to a standing 
accretion shock, and the activities of convective overturns and SASI become vigorous with 
time behind the standing shock. Since anisotropies of the neutrino flux and matter motions 
in the postbounce phase are governed by the non-linear hydrodynamics, GWs emitted in 
this epoch also change stochastically with time (indicated by "Stochastic" in the table, e.g., 
section [3TT1 for more details). For detecting these signals for a Galactic CCSN event with a 
good signal-to-noise ratio, we need next generation detectors such as the advanced LIGO, 
LCGT, and the Fabry- Perot (FP)-type DECIGO. Neutrino signals at the non-linear epoch 
also change stochastically with time , which is expected to be detected by currently running 
detector such as by IceCube and SK (IBrandt et al.ll201ll ; iMarek et al.ll2009t iLund et al.ll201Q[ ) 
for the galactic event. It should be also mentioned that detailed neutrino transport includ- 
ing the most up-to-date interaction rates should be accurately implement ed in the multi- 
D sim u lations to make a r eliable prediction of the neutrino signals (e.g., iHudepohl et al. 



(120101 ); Fischer et al.l (120101 ) ). When material behind the stalled shock successfully absorbs 
enough neutrino energy to be gravitationally unbound from the iron core, the up process is 
expect ed to successfully explain some light proton- rich nuclei in the neutrino-driven wind 
phase (IPruet et al.l 120061 ; iFrohlich et al.l 120061 ; IWanajo fc Jankal 1201 ll ) . When the revived 
supernova shock (successfully) attains the kinetic energy as energetic as 10 51 erg, explosive 
nucleosynthesis that proceeds in the SASI-aided low-mode explosions is expected to partly 
explain the solar abundance yields as well as the observed non-uniform morphology of the 
synthesized elements (e.g., section I2".3.ip . 

One important notice here is that explosion energies obtained in the state-of-the-art 
multi-D simulations are typically underpowered by one or two orders-of-magnitudes to ex- 
plain the canonical supern ova kinetic energy ( y 10 51 erg, e.g., Table 1). Moreover, the 



softer EOS, such as of the lLattimer fc Swestyl (Il99ll ) (LS) EOS with an incompressibility 
at nuclear densities, K, of 180 MeV, is employed in those simulations (e.g., in multi-D sim- 
ulations of the MPA, Oak Ridge+, and Tokyo-I- groups in Table 1). On top of a striking 
evidence that favors a stiffer EOS based on the nuclear experimental data (K = 240 ± 20 
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neutron star of ~ 2M^ ( 



in 



O'Connor & Ott 2011 



MeV, IShlomo et al.l (120061) ). the soft EOS may not account for the recently observed massive 



Demorest et al. 2010) (see the maximum mass for the LS180 EOS 



Kiuchi fc Kotakd 120081). With a stiffer EOS, the explosion energy 



may be even lower as inferred from lMarek fc Jankal ( 20091 ) who did not obtain the neutrino- 
driven explosion for their model with K = 263 MeVI 32 !. What is then missing furthermore ? 
We may get the answer by going to 3D simulations (section 2.1) or by taki ng into account new 



ingredients, such a s exotic physics in the core of the protoneutron star (ISagert et al.l 12009 



Fischer et al.l 1201 lh . viscous heating by the magneto-rotat ional instability (section |3.0.3|) 
or energy dissipation via Alfven waves (ISuzuki et al.ll2008l ). In additio n to the 3D effects 



GR is expected to h elp the onset of multi-D neutrino-driven explosions (IMueller et al.l 12011 



Kuroda et al.ll2012f l aa l At present, the most up- to-date neutrino transport code can treat the 
multi-energy and multi-angl e transport in 2D (jOtt et al.ll2008l ) and even in 3D simulations 
(jSumiyoshi &: Yamadall2012l ) but only limited to the Newtonian simulations. Extrapolating 
a current growth rate of supercomputing power, the dream simulation (i.e. 6D simulations 
with full Boltzmann neutrino transport in full GR simulations) is likely to be practicable 
not in the distant future (in our life-time!) presumably by using the exascale platforms. 

In addition to these numerical advancements, the physical understanding of the su- 
pernova theory via analytical studies is also in a steady progress. What determines the 
saturation levels of SASI ? A careful analysis on the parasitic instabilities has been reported 
to answer this question (IGuilet et al.l l2010l ). What determines mode couplings between 
small- scale convection eddies a nd large-scale SASI modes ? To apply the theory of turbu- 
lence (IMurphy &: Meakinll201ll ) should put a milestone to address this question. Besides the 
two repre sentative EOSs of Lattimer-Swesty and Shen, new sets of EOSs have been recently 
reported (IHempel fc Schaffner-Bielichll2010l ; iFurusawa et al.l 1201 ll : IShen et al.l 1201 11 ). Some 
of them will enable SN modellers to go beyond a single nucleus approximation, which is an 
important improvement for an accurate description of neutrino-nucleus interaction. 

Rotation, albeit depending on its strength, can give a special direction (i.e., spin axis) 
in the supernova cores. As discussed in section |2.2.2[ stochastic nature of the GWs becomes 
weak in the presence of rotation as a result of the spiral SAsjj^|( indicated by "excess for 



32 0n the other hand, they obtained 2D explosions for Shen EOS (K 
communication) . 

33, 



281MeV) (H-T. Janka, private 



J e.g. 



Mueller et al.l ([20121 ) in 2D simulations but with detailed neutrino transport or lKuroda et al.l (|2012l ) 



in 3D simulations but with a very approximate neutrino transport 

34 It should be noted that the fea ture can be seen even in a slowly rotating progenitor that agrees with 
recent stellar evolution calculation ( Maeder fc Mevnet 2000t Heger et al. 2005 ) (see section 12.2.21 for more 
details). 
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equator" in Table 2). Although it is not easy to detect these low-frequency GWs from 
anisotropic neutrino emission by currently running laser interferometers, a recently proposed 
space-based interferometers (like FP-type DECIGO) would permit detection for a galactic 
event. Contributed by the neutrino GWs in the lower frequency domain, the total GW 
spectrum is expected to become rather flat over a broad frequency range below ~ 1 kHz. 
This GW feature obtained in the context of the SASI-aided neutrino-driven mechanism is 
different from the one in the other candidate super nova mechanisms , such as the MHD 
mechanism (section [XT]) and the acoustic mechanism (jOtt et al.ll2006al ). thus could provide 
an important probe into the explosion mechanism. Rapid rotat ion induces a polar excess of 



neutrino emission w i th a harder spectrum than on the equator (IJanka fc Monchmeyerlll989 



Kotake et al.l l2003at iBrandt et al.ll201ll ) (indicated by "polar excess" in the table). This 



should be also the case of the MHD mechanism, which deserves further investigation for 
a quantitative discussion (symbolized by Polar excess (?) in the table). These (rotation- 
induced) neutrino signatures, if observed, will carry an important information about the 
angular momentum profiles hidden deep inside massive stellar cores. Concerning the z/p 
process and explosive synthesis in rapidly rotating cores (as well as in the MHD explosions), 
there still remains a vast virgin territory for further investigation (symbolized by " z/p process 
(?)" in the table). 

If the neutrino-heating mechanism (or some other mechanisms) fails to blow up m assive 
stars, central PNSs collapse to BHs. Recent GR simulations by lOtt et al.l ( 1201 ll ) show 
that the significant GW emission is associated at the moment of the BH formation, which 
ca n be a promising targe t of the advanced LIGO for a galactic source. As pointed out 
by ISumiyoshi et al.l (120061 ) , the disappearing neutrino signals that mark the epoch of BH 
formations also can be a target of SK and IceCube (llceCube Collaboration et al.l l201ll ). 
When quarks and hyperons appear in the postbounce core, a neutrino burst pro duced by 
the sudden EOS softening and by the subseque nt rebounce dSaeert et al.l 120091) is likely 
to be detected by IceCube for a galactic event (llceCube Collaboration et al.l 1201 ll ) . The 
interval between core-bounce and the BH formation depends on the details of exotic physics 
in the super-dense PNS cores iNakazato et aliboiOal lbL All of these observational signatures 
should provide an important probe in t o the so-called dense QCD region in the QCD phase 



diagram (e.g.. iFukushima &: Hatsudal (120111 ) for recent review) to which lattice calculations 



are hardly accessible at present. Concerning photons, no optical outbursts are expected in 
the BH-forming SNe (indicated by "no photon (?)" in Table 2), if not for rapid rotation and 
strong magnetic fields prior to core-collapse. 

If a precollapse core rotates enough rapidly (typically initial rotation period less than 
4 s) with strong magnetic fields (higher than ~ 10 11 G), the MHD mechanism can produce 
bipolar explosions along the rotational axis predominantly driven by the field wrapping 
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processes (section |3.0.2[) . If the MRI can be sufficiently resolved in global simulations, the 
MRI would exponentially amplify the initial magnetic fields to a dynamically important 
strength within several rotational periods (section 13.0.31 and I3.0.4p . The next-generation 
supercomputing resources are needed again to see the outcome. If such simulations would be 
executable, the MHD outflows will be produced even for more weakly magnetized cores than 
currently predicted. The GW signals in the MHD explosions are characterized by a burst-like 
bounce signal plus a secularly growing tail (section I3.1[ indicated by "burst and tail" in the 
table). Similar to the neutrino GWs in the neutrino-driven mechanism, a future detector 
(like FP-type DECIGO) is again necessary for detecting the low-frequency tail component. 
As discussed in section 13.2.11 the pole to equator anisotropy of the shock propagation in the 
MHD explosions affect the neutrino signals through the MSW effects. The anisotropic shock 
passage to the H-resonance regions leads to a sudden decrease in the SK events (~ 2500 
events for a galactic source), which is quite possible to identify by the SK-class detectors 
(indicated by "anisotropy in SK events" in Table 2). Since the MHD mechanism has such 
distinct signatures, a planned joint analysis of neutrino and GW data ( iLeonor et al.l 120101 ) 
would be potentially very powerful to tell the MHD mechanism from the other candidate 
mechanisms. If the resonant spin-flavor conversion (indicated as RSF in Table 2) occurs in 
the highly magnetized core, the neutronization burst of v e converts to t hat of whic h is th us 



expected to be a probe into the magnetic moment of neutrinos (e.g., lAndo fc S ato (2003)). 



Finally the MHD explosion s are considered to be a possible r-process cite (INishimura et al. 



20061 ; iFujimoto et al.ll200 7l 35 l This is because the mass ejection from the iron core occurs in 
much shorter timescales compared to the delayed neutrino-heating mechanism, so that the 
ejecta can stay neutron rich before it becomes proton-rich via neutrino capture reactions. 
Most of these possibilities regarding the MHD mechanism have been proposed so far in 
simulations with a crude treatment of neutrino transport, which is now rather easily re- 
examined by the state-of-the-art multi-D simulations (if the code is MHD). 

Finally, what is about the story if the MHD mechanism fails ? (indicated by "path 
to hypernovae in Table 2). The rapidly rotating PNS collapses into a BH, probably lead- 
ing to the formation of accretion disk around the BH. Neutrinos emitted from the accre- 
tion disk heat matter in the polar funnel region to launch outflows or strong magnetic 
fields in the cores of the order of 10 15 G also play an active role both in driving the 
magneto-drive n jets a nd in extracting a sign i ficant amount of energy f rom t he BH (e.g., 
Wheeler et aD (boooh ; iThompson et al.l ( 120041 ); llJzdenskv fc MacFadven! (l2007h and see ref- 
erence s therein). This p i cture , often referred as collapsar (e.g., iMacFadyen fc Woosley 



(119991 ); IMacFadyen et al.l ( 120011 )). has been the working hypothesis as a central engine of 



See collective references in lWanaio fc Jankal (| 2 1 lh for other plausible r-process cites 
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other candidate mechanism including magnetar models). However it is still controversial 
whether the generation of the relativistic outflows predominantly proceeds via the neutrino- 
heating mechanism or the MHD mechanism. In contrast to a number of findings illustrated 
in Table 2, much little things are known about the BH- forming supernovae and also about the 
collapsar. This is mainly because the requirement for making a realistic numerical modeling 
is very computationally expensive, which (at least) necessitates the multi-D MHD simula- 
tions not only with GR for handling the BH formation (thus full GR), but also with the 
multi-angle neutrino transfer for treating highly anisotropic neutrino radiation from the ac- 
cretion disk. To get a unified picture of massive stellar death, we need to draw a schematic 
picture (like Table 2) also in the case of the BH-forming supernovae. We anticipate that our 
long-lasting focus (and experience) on the explosion dynamics of canonical CCSNe will be 
readily applicable to clarifying the origin of the gigantic explosion energy of hypernovae and 
also unraveling (ultimately) the central engine of long-duration GRBs. This should provide 
yet another grand challenge in computational astrophysics in the next decades. 

As repeatedly mentioned so far, all the numerical results in this article should be tested 
by the next-generation calculations by which more sophistication is made not only in the 
treatment of multi-D radiation transport (both of neutrino and photon), but also in multi-D 
hydrodynamics including stellar rotation and magnetic fields in full GR. From an optimistic 
point of view, our understanding on every issue raised in this article can progress in a step- 
by-step manner at the same pace as our available computational resources will be growing 
bigger and bigger from now on. Since 2009, several neutrino detectors form the Super- 
nova Early Warning Systems (SNEWS) to broadcas t the alert to astronomers to let them 
know the arrival of neutrinos (lAntonioli et al.l 120041 ) . Currently, Super-Kamiokande, LVD, 
Borexino, and IceCube contributes to the SNEWS, with a number of other neutrino and 
GW detectors planning to join in the near future. This is a very encouraging news towards 
the high-precision multi-messenger astronomy. The interplay between the detailed numeri- 
cal modeling, the advancing supercomputing resources, and the multi-messenger astronomy, 
will remain a central issue for advancing our understanding of the theory of massive stellar 
core-collapse for the future. A documentary film recording our endeavours to make practica- 
ble the "multimessenger astronomy of CCSNe" seems not to show "fin" immediately and is 
becoming even longer as far as the three components evolve with time. To raise the edifice, 
it is becoming increasingly important to bring forward a world-wide, multi-disciplinary col- 
laboration among different research groups. We believe that such an approach would provide 
the shortest cut to get a deeper understanding of a number of unsettled and exciting issues 
that we were only able to touch in this article. 
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